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1.  INTRODUCTION 


Assignments  for  aeronaut, cal  radio  in  the.  radio  frequency  spectrum  must  be  made 
so  as  to  provide  reliable  services  for  an  increasing  air  traffic  density  [17]"*. 
Potential  interference  between  facilities  operating  on  the  same  or  on  adjacent  chan¬ 
nels  must  be  considered  in  expanding  present  services  to  meet  future  demands.  Ser¬ 
vice  quality  depends  on  many  factors,  including  signal  strength  and  the  desired-to- 
undesired  signal  ratio  at  the  receiver.  These  parameters  vary  with  receiver  location 
and  time  even  when  other  parameters,  such  as  antenna  gain  and  radiated  powers,  are 
fixed. 

In  1973,  an  air/ground  propagation  model  developed  at  the  Department  of  Commerce 
Boulder  Laboratories  (DOC-BL)  by  the  Institute  for  Telecommunication  Sciences  (ITS) 
for  the  Federal  Aviation  Administration  (FAA)  was  documented  in  detail  [14].  This 
IF-73  ( ITS-FAA- 1973)  propagation  model  has  evolved  into  the  IF -77  ( ITS-FAA-1 977 ) 
model,  which  is  applicable  to  air/air,  air/ground,  air/satellite,  ground/ground,  and 
ground/satellite  paths.  The  If -77  has  been  incorporated  into  a  number  of  computer 
programs  that  are  useful  in  estimating  the  service  coverage  of  radio  systems  oper¬ 
ating  in  the  frequency  band  from  0.1  to  20  GHz.  These  programs  may  be  used  to 
obtain  a  wide  variety  of  computer-generated  microfilm  plots  [19,20].  Extensive 
comparisons  of  IF-77  predictions  with  measured  data  have  been  made  [21],  and  an 
atlas  of  basic  transmission  loss  predictions  was  generated  using  the  model  [8,22], 

The  previously  published  documentation  for  IF-77  covered  only  the  extensions 
[15]  made  to  IF-73  and  relied  on  IF-73  documentation  [14]  to  cover  most  of  the 
model.  This  approach  was  used  to  highlight  the  extensions  and  simplify  the  documen¬ 
tation  process.  However,  an  integrated  description  of  the  model  is  needed  to  give 
potential  users  a  better  understanding  of  the  overall  model.  Such  a  description  is 
provided  in  this  report.  Although  most  of  the  previous  descriptions  are  repeated 
here  in  an  integrated  form,  some  items  are  simply  referenced.  Details  concerning 
computer  programs  are  not  included,  but  it  is  anticipated  that  suCh  detail  for  a 
few  specific  programs  will  be  made  available  in  later  reports. 

Except  where  otherwise  indicated,  all  equations  provided  here  are  dimensionally 
consistent;  e.g.,  all  lengths  in  a  particular  equation  are  expressed  in  the  same 
units.  Distances  and  heights  are  always  in  kilometers,  angles  always  are  in  radians, 
and  frequency  is  always  in  megahertz.  In  Section  2.3,  wavelength  is  in  meters,  else- 


'Referencos  arc  listed  alphabetical iy  by  author  at  the  end  of  the  report  so  that 
reference  numbers  do  not  appear  sequentially  in  the  text. 


where  it  is  in  kilometers.  Braces  are  used  around  parameter  dimensions  when  partic¬ 
ular  units  are  called  for  or  when  a  potential  dimension  difficulty  exists.  A  list, 
of  symbols  is  provided  in  the  appendix. 

2.  PROPAGATION  MODEL 

The  IE-77  propagation  model  is  applicable  to  air/ground,  air/air,  ground/satel¬ 
lite,  and  air/satellite  paths.  It  can  also  be  used  for  ground/ground  paths  that  are 
1 ine-of-sight  or  smooth  earth.  Model  applications  are  restricted  to  telecommunica¬ 
tion  systems  operating  at  radio  frequencies  from  about  0.1  to  20  GHz  with  antenna 
heights  greater  than  0.5  m.  In  addition,  radio-horizon  elevations  must  be  less  than 
the  elevation  of  the  higher  antenna.  The  radio  horizon  for  the  higher  antenna  is 
taken  either  as  a  common  horizon  with  the  lower  antenna  or  as  a  smooth  earth  horizon 
with  the  same  elevation  as  the  lower  antenna  effective  reflecting  plane  (Sec.  3.3). 

At  0.1  to  20  GHz,  propagation  of  radio  energy  is  affected  by  the  lower,  nonion- 
ized  atmosphere  (troposphere),  specifically  by  variations  in  the  refractive  index 
of  the  atmosphere  [1,  2,  3,  4,  6,  9,  18,  23,  30,  33,  34],  Atmospheric  absorption 
and  attenuation  or  scattering  due  to  rain  become  important  at  SHF  [23,  Ch.  8;  28; 

34,  Ch.  3].  The  terrain  along  and  in  the  vicinity  of  the  great  circle  path  between 
transmitter  and  receiver  also  plays  an  important  part.  In  this  frequency  range, 
time  and  space  variations  of  received  signal  and  interference  ratios  lend  themselves 
readily  to  statistical  description  [18,  25,  27,  32,  34,  Sec.  10]. 

Conceptually,  the  model  is  very  similar  to  the  Longley-Rice  [26]  propagation 
model  for  propagation  over  irregular  terrain,  particularly  in  that  attenuation 
versus  distance  curves  calculated  for  the  1 ine-of-sight  (Sec.  5),  diffraction 
(Sec.  4),  and  scatter  (Sec..  6)  regions  are  blended  together  to  obtain  values  in 
transition  regions.  In  addition,  the  Longley-Rice  relationships  involving  the 
terrain  parameter  Ah  are  used  to  estimate  radio-horizon  parameters  wnen  such  infor¬ 
mation  is  not  available  from  facility  siting  data  (Sec.  3.2).  The  model  includes 
allowance  for: 

1.  Average  ray  bending  (Sec.  3). 

2.  Horizon  effects  (Sec.  3). 

3.  Long-term  power  fading  (Sec.  10.1). 

4.  Antenna  pattern  (elevation  only)  at  each  terminal  (Sec.  5.2.4). 

5.  Surface  reflection  multipath  (Sec.  10.2). 
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6.  Tropospheric  multipath  (Sec.  10.3). 

7.  Atmospheric  absorption  (Sec.  9). 

8.  Ionospheric  scintillations  (Sec.  10.5). 

9.  Rain  attenuation  (Sec.  10.4). 

10.  Sea  state  (Sec.  5.2 .2). 

11.  A  divergence  factor  (Sec.  5.2.1). 

12.  Very  high  antennas  (Sec.  8). 

13.  Antenna  tracking  options  (Sec.  5.2.4). 

Computer  programs  that  utilize  IF-77  calculate  transmission  loss,  power  avail¬ 
able,  power  density,  and/or  a  desired-to-undesired  signal  ratio.  These  parameters 
are  discussed  in  Sections  2.1  through  2.4.  A  computational  flow  chart  is  provided 
in  Section  2.1  for  transmission  loss. 

2.1  Transmission  Loss 

Transmission  loss  has  been  defined  as  the  ratio  (usually  expressed  in  decibels) 
of  power  radiated  to  the  power  that  would  be  available  at  the  receiving  antenna 
terminals  if  there  were  no  circuit  losses  other  than  those  associated  with  the 
radiation  resistance  of  the  receiving  antenna  [34,  Sec.  2],  Transmission  loss 
levels,  L(q),  that  are  not  exceeded  during  a  fraction  of  the  time  q  (or  100  q  percent 
of  the  time)  are  calculated  from: 

L(q)  =  Lb(0.5)  +  Lgp  -  -  GER  -  Y£(q)  dB  (1) 

where  Lb(0.5)  is  the  median  basic  transmission  loss  [32,  Sec.  10.4],  Lgp  is  the  path 
antenna  gain-loss  [26,  Sec.  1-3],  G^y  and  G^R  are  free-space  antenna  gains  for  the 
transmitter  and  receiver  at  the  appropriate  elevation  angle,  respectively,  and  Y„(q) 
is  the  total  variability  from  equation  229  of  Section  10;  i.e.:  (229). 

Median  basic  transmission  loss,  Lb(0.5),  is  calculated  from: 

Lb(°.5)  -  Lbr  +  Ay  +  Aa  dB  (2) 

where  Lbr  is  a  calculated  reference  level  of  basic  transmission  loss,  Ay,  is  a  con¬ 
ditional  adjustment  factor,  and  A,  is  atmospheric  absorption  from  (228)  of  Section  9. 
The  factor.  Ay,  from  (243)  jf  Section  10.1,  is  used  to  prevent  available  signal 


power;;  from  exceeding  levels  expected  for  free-space  propagation  by  an  unrealistic 
amount  when  the  variability  about  1.^(0.  5)  is  large  and  Lb(0.5)  is  near  its 
free-space  level. 

Free-space  basic  transmission  loss,  L^,  from  (226)  of  Section  8,  terrain  atten¬ 
uation,  Ay,  from  (221)  of  Section  7,  and  a  variability  adjustment  term,  Ve(0.5), 
from  (240)  of  Section  10.1,  are  used  to  determine  Lbr;  i.e.: 

Lbr  =  Lbf  +  AT  -  Ve(0‘5)  dB  0) 

The  value  for  L  in  (1)  is  taken  as  0  dB  in  the  IF-77  model.  This  is  valid 
when  (a)  transmitting  and  receiving  antennas  have  the  same  polarization  and  (b)  the 
maximum  antenna  gain  is  less  than  50  dB  [26,  Sec.  1-3]. 

Values  of  G^y  and  GyR  for  (1)  are  obtained  from  G^y  =  Gy  +  G^y  and 
GER  =  GR  +  GNR;  i,e* : 


GET,R  =  GT,R  +  GNT,R  dBl  '4) 

where  Gy  R  is  the  transmitter  or  receiver  main-beam  maximum  free  space  antenna  gain 
in  decibels  greater  than  isotropic  (dBi),  and  GNy  R  is  a  normalized  transmitting  or 
receiving  antenna  gain  in  decibels  greater  than  maximum  gain  that  gives  relative 
gain  for  the  appropriate  elevation  angle.  These  gains  are  all  model  input  parameters 
Normalized  vertical  (elevation)  antenna  patterns  are  used  to  define  G^y  R  so  that 
gain  values  can  be  obtained  for  elevation  angles  at  which  the  maximum  gain  is  not 
appropriate.  The  calculation  of  these  elevation  angles  is  discussed  in 
Section  5.2.4.  Horizontal  antenna  patterns  are  :<t  usually  considered  part  of  the 
IF-77  model,  but  an  allowance  for  them  can  be  made  by  adjusting  Gy  R  values. 

However,  horizontal  patterns  have  been  included  in  one  computer  program  called  TWIRL 

[20,  p.  61]. 

Variabilities  associated  with  long-term  power  fading,  surface  reflection  multi- 
path,  tropospheric  multipath,  rain  attenuation,  and  ionospheric  scintillation  are 
included  in  Y^(q)  of  (1).  These  are  all  discussed  in  Section  10.  Since  the  adjust¬ 
ment  terms  associated  with  long-term  power  fading  are  included  in  the  calculation 
of  Lb(0.5),  YJC(0.  5)  =  0. 

A  computational  flow  diagram  for  L(q)  is  provided  in  Figure  1.  Although  only 
those  equations  discussed  in  this  section  are  included  in  the  diagram,  section  ref¬ 
erences  are  provided  for  the  other  calculations  involved.  Horizon  and  diffraction 


FIGURE  1.  COMPUTATIONAL  FLOW  DIAGRAM  FOR  L(q). 


5 


region  parameters  are  determined  prior  to  calculation  of  since  these  parameters 
are  used  to  define  region  boundaries. 

2.2  Power  Available 

Power  available  as  calculated  in  IF-77  is  taken  as  the  power  available  from  the 
receiving  antenna  terminals  under  matched  conditions  when  internal  heat  losses  of 
the  receiving  antenna  and  path  antenna  gain- loss  are  neglected.  Compensation  for 
internal  heat  loss  or  gain-loss  factors  needed  to  refer  the  available  power  to  some 


point  in  the  receiving  system  other  than  the  receiving  antenna  terminals  can  be  made 
by  an  appropriate  adjustment  to  the  radiated  power  or  antenna  gains  used  for  computer 
program  input. 

Power  available  P  (q)  levels  exceeded  for  a  fraction  of 
using  Lb(0.5)  from  (2),  G^j  R  from  (4),  and  Y>.(q)  from  (229) 

time  q  are  determined 

of  Section  10  with: 

Pa(q)  =  EIRPG  +  Gnt  +  Gnr  -  Lb(0.5)  +  Y>:(q)  dBW 

(5) 

EIRPG  =  EIRP  +  Gr  dBW 

(6) 

EIRP  =  PTR  +  Gt  dBW 

(7) 

Here  EIRP  is  equivalent  isotropically  radiated  power;  PTR  in  decibels  greater  than 
1  W  (dBW)  is  the  , total  power  radiated  by  the  transmitting  antenna;  and  Gj  R  in  deci¬ 
bels  greater  than,  isotropic  (dB i )  is  the  maximum  gain  of  the  transmitting  antenna  or 
receiving  antenna,  respectively  (Sec.  2.1).  Losses  (e.g.,  line  losses)  associ¬ 
ated  with  the  transmitting  system  should  be  considered  in  calculating  radiated 
power  from  transmitter  output  power.  Normal. zed  antenna  gains  ( and  G^R)  in 
decibels  greater  than  maximum  gain  (G^  or  GR)  are  included  in  (5)  to  allow  for 
antenna  directivity  when  maximum  gains  are  not  appropriate  (i.e.f  the  cmennas 
are  not  pointed  at  each  other).  Computer  programs  utilizing  IF-77  have  tracking 
options  that  allow  antennas  to  track  each  other  (Sec.  5.2.4). 


2.3  Power  Density 

Power  density  SR(q)  in  decibels  greater  than  1  watt  per  square  meter  that  is 
exceeded  for  a  fraction  cf  the  time  q  is  determined  using  the  parameters  discussed 
previously  along  with  the  effective  area,  Aj .  of  an  isotropic  antenna;  i.e,: 

SR(q)  =  EIRP  +  Gnt  -  Lb(0.S)  +  Y?.(q)  -  Aj  dR-W/sq  m  (8) 
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(9) 


Values  of  Aj  are  determined  from: 

Aj  =  10  log  (A^/4ir)  dB-sq  m 
where  Am  is  the  wavelength  in  meters  (33,  Sec.  4.11).  For  a  frequency  of  f  [MHz]: 

Am  =  299. 7925/f  m  (10a) 

A  =  A  71000  km  ( 1  Ob ) 

m  ' 


2.4  Desired-to-Undesired  Signal  Ratio 

Desired- to-undesired  signal  ratios  that  are  available  for  at  least  a  fraction 
time  q,  D/U(q)  dB,  at  the  terminals  of  a  lossless  receiving  antenna  are  calculated 
using  [11,  Sec.  3]: 

D/U(q)  *  D/U(0.5)  +  YDU(q)  dB  (11) 

The  median  value  of  D/U(0. 5)  and  the  variability  YRU(q)  of  D/U  are  calculated  as: 
D/U(0. 5)  -  [Pa(0.5)]Desjred  -  CPa(0.5)]Undes1red  02) 

Wq^  *  -  y  Oj:(q)-l2Uesjred  +  [VE(1  -  <l)32ur,des1red  dB  (13) 

-  for  q  >  0.5 
+  otherwise 

Values  of  P  (0.5)  are  calculated  from  (5)  where  Y^(0.5)  =  0  by  using  parameters 
appropriate  for  either  the  desired  or  undesired  facility.  Applicable  variabilities 
are  calculated  using  the  methods  described  in  Section  10.  Note  that  Y^y(q)  requires 
the  undesired  facility  for  (1  -  q);  e.g.: 

YDU ^ °‘ S)  =  "  y  [ Y>; ( o. B )] ‘''Desired  +  2^2|Jndesired 
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3.  HORIZON  GEOMETRY 

Calculations  associated  with  horizon  geometry  involve  the  use  of  the  effective 
earth  radius  concept  in  which  ray  bending  caused  by  refraction  within  the  troposphere 
is  simplified  by  using  straight  rays  above  an  earth  with  an  effective  radius  that  is 
selected  to  compensate  for  the  ray  bending  [3,  Sec.  3.6;  30].  The  effective  earth 
earth  radius,  a,  is  calculated  [34,  Sec.  4]  using  the  minimum  monthly  mean  surface 
refractivity  referred  to  mean  sea  level,  N  ,  and  the  height  of  the  effective  reflec¬ 
tion  surface  above  mean  sea  level  (msl),  hr:  i . e . : 


N$  =  greater  of 


NQ  exp  (-0.1057  h  ) 
or 
?00 


a^  =  6370  km 
o 

a  =  a  [1 '  -  0.04665  exp  (0.005577  N^.)]-1 


M-uni ts 


km 


(14) 

(15) 

(16) 


Here  Ng  is  surface  refractivity  at  the  effective  reflecting  surface,  and  aQ  is  the 
actual  earth  radius  to  three  significant  figures.  Both  NQ  [20,  p.  74,  p.  94]  and 
hr  [20,  p.  83]  are  model  input  parameters. 

When  high  (>T  km)  antennas  are  involved,  geometry  based  solely  on  the  effective 
radius  method  may 'overestimate  ray  bending  so  that  smooth  earth  horizon  distances 
become  excessive  [31].  This  difficulty  is  compensated  for  in  IF-77  by  the  use  of 
ray  tracing  in  the  determination  of  some  key  parameters  such  as  effective  terminal 
altitudes  (Sec.  3.1),  smooth  earth  horizon  distances  (Sec.  3.1),  and  effective  dis¬ 
tance  (Sec.  10.1).  In  IF-77,  ray  tracing  is  performed  through  an  exponential 
atmosphere  [3,  equations  3.43,  3.44,  3.40]  in  which  the  refractivity,  N,  varies  with 
height  above  msl,  h,  as: 

N  -  N$  exp  [-  C  (h  -  h(J]  N-units  (17) 

where 


and 

AN  -  -7.32  exp(0. 006677  N  )  N-units/km  (19) 
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Thayer's  algorithm  [37]  for  ray  tracing  through  a  horizontally  stratified  atmosphere 
is  used  with  layer  heights  (above  hr)  taken  as  0.01,  0.02,  0.05,  0.1,  0.2,  0.305, 

0.5,  0.7,  1,  1.524,  2,  3.048,  5,  7,  10,  20,  30.48,  50,  70,  90, .110,  225,  350,  and 
475  km.  Above  475  km  ray  bending  is  neglected;  i.e.,  rays  are  assumed  to  be  straight 
relative  to  a  true  earth  radius,  a0.  The  computer  subroutine  used  for  ray  tracing 
[14,  Sec.  B.4.1,  RAYTRAC]  was  written  so  that:  (a)  the  initial  ray  elevation  angle 
may  be  negative;  (b)  if  the  initial  angle  is  too  negative,  it  will  be  set  to  a  value 
that  corresponds  to  grazing  for  a  smooth  earth;  and  (c)  the  antenna  heights  may  be 
very  large,  e.g.,  satellites. 

The  IF-73  model  was  actually  developed  for  transmission  from  a  ground-based 
facility  to  an  airborne  receiver;  but,  because  of  reciprocity,  it  could  also  be  used 
for  air-to-ground  transmissions.  Although  the  IF-77  model  does  not  require  that  a 
terminal  be  ground-based,  or  that  a  terminal  be  airborne,  we  will  nevertheless  refer 
to  the  lowest  terminal  as  the  facility  and  the  higher  terminal  as  the  aircraft. 
Furthermore,  the  facility  is  taken  as  terminal  1  so  that  variables  associated  with 
it  have  a  1  in  their  subscripts;  and  the  aircraft  is  taken  as  terminal  2. 

Horizons  for  both  terminals  are  determined  by  the  input  parameters  used  to  de¬ 
fine  the  facility  horizon;  i.e.,  the  aircraft  horizon  is  either  taken  as  the  facility 
horizon  obstacle  (common  horizon)  or  as  a  smooth  earth  when  the  aircraft's  view  of 
the  facility  obstacle  would  be  shadowed  by  a  smooth  earth.  The  calculation  of  air¬ 
craft  horizon  parameters  (Sec.  3.3)  involves  prior  determination  of  effective 
antenna  heights  (Sec.  3.1)  and  the  facility  horizon  (Sec.  3.2). 

3.1  Smooth-Earth  horizons 

Antenna  height  information  is  input  to  IF-77  via  the  variables  hr,  hsl ,  hc, 
h$mp  and  h^  where  these  variables  are  defined  in  Figure  2  as  in  the  Application 
Guide  [20,  pp.  80,  81,  83,  88;  101].  Antenna  heights,  above  the  effective 

reflection  surface  are  obtained  from  these  by  using: 

hrl  =  hsl  +  hsml  -  hr  k,n 
hr2  =  h2  "  hr  km 

Antenna  heights  above  mean  sea  level  (ins  1 }  are  given  by: 

hl  ,2  =  hrl ,2  +  hr  km  ^21 ^ 


(20a) 

(20b) 
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AIRCRAFT  ALTITUDE  ABOVE  msl 


FACILITY  ANTENNA  HEIGHT  ABOVE  fss 


FACILITY  ANTENNA  COUNTERPOISE  ABOVE  fss 


FACILITY  SITE  SURFACE  (fss)  ELEVATION  ABOVE  msl 


EFFECTIVE  REFLECTION  SURFACE  ELEVATION  ABOVE  msl 


MEAN  SEA  LEVEL  (msl) 


VALID  INPUT  CONSTRAINTS 
0<  hr  £  4  km 
°ihsm,  i  4  km 


0.0005  km  < 

h,  <  h2 


Note  that  aircraft  altitude  is  elevation  above  msl  while 
the  facility  antenna  height  is  measured  with  respect  to  fss. 


FIGURE  2.  INPUT  ANTENNA  HEIGHTS  AND  SURFACE  ELEVATIONS  FOR  IF-77. 
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The  height  of  the  facility  antenna  above  the  counterpoise  is  given  by: 


hfc  “  '’si  “  hc  km  ^22^ 

When  the  effective  earth  geometry  (straight  rays)  overestimates  bending,  effective 
antenna  heights,  hg^  2«  are  taken  as  heights  that  will  yield  the  smooth-earth  hori¬ 
zon  distances  obtained  via  ray  tracing  when  used  with  effective  earth  geometry. 

As  illustrated  in  Figure  3,  effective  heights  are  lower  than  actual  heights. 


FIGURE  3.  EFFECTIVE  HEIGHT  GEOMETRY . 


The  IF -77  model  uses  ray  tracing  to  determine  smooth-earth  horizon  ray  dis¬ 
tances  associated  with  both  terminals,  dLsR1  2,  that  are  used  in  the  calculation  of 
effective  antenna  heights.  These  distances  are  determined  by  ray  tracing  from  the 
effective  reflecting  surface  elevation  of  the  earth's  surface  to  the  respective 
antenna  heights.  The  initial  takeoff  angle  used  is  0"  and  the  surface  ref racti vi ty , 
Nc,  is  calculated  from  NQ  with  (14).  Values  for  dL$R1  2  and  a  from  (16)  are  used  to 
determine  the  difference  in  actual  and  effective  antenna  heights.  Ah  ,  as  follows 


km  (24) 


h  ,  o  "  lesser  of 
el  ,2 

°'5  dLsRl,2/a  if  9sR1 ,2  -  0,1  rad 

m 

a[sec(e$R1  2)  -  1]  otherwise 

Ahel,2  =  hrl ,2  "  hel,2  km 


(25) 


The  final  value  of  a  smooth  earth  horizon  distance,  d^  2>  is  taken  as  the 
ray  tracing  value  for  high  antennas  (Ah^  2  >  0)  or  computed  via  effective  earth 
radius  geometry,  dLsE1  2;  i.e. : 


dLsEl  ,2 

= 

V  *el , 

2  km 

(26) 

dURl,2  if 

Ahel2>0 

dLsl  ,2 

= 

km 

(27) 

dLsEl  ,2 

otherwise 

dLs 

=  d 

Lsl  +  dLs2 

km 

(28) 

In  addition,  the  ray  elevation  angle  Qes^  ^  resulting  from  ray  tracing  is  taken  as 
the  final  ray  elevation  angle  when  Ah^  2  >  0;  i.e.: 


si, 2 

esl  ,2 


dLsl ,2/a  rad 


0esRl ,2  lf  Ahel,2  >  0 


otherwise 


rad 


(29) 


(30) 


3.2  Facility  Horizon 


The  IF-73  model  allowed  the  facility  horizon  to  be  specified  by:  (a)  any  two 
horizon  parameters  (elevation,  elevation  angle,  or  distance);  (b)  estimated  with 
any  one  horizon  parameter  and  the  terrain  parameter.  Ah;  (c)  estimated  from  Ah 
alone;  or  (d)  calculated  for  smooth-earth  conditions  [14,  Fig.  14].  Figure  4  illus¬ 
trates  the  facility/ horizon  geometry  involved  in  the  IF-73  formulation.  Horizon 
distances,  d^ ;  horizon  elevation,  h^  ,  and  horizon  elevation  angle,  0  ^ ;  Tor 
effective  earth  geometry  are  related  to  each  other  by: 


0  -  Tnn  -1  1  hlE1  "  h]  d|“E1  1  nd 

9eEl  ‘  Tan  (  dLE,  '  ~W )  rad 

(31) 

(dLE,)2 

hLEI  '  hl  +  2a  +  dLEl  tan  0eEl  km 

(32) 

dLEl  =  -^2a^hLEl  ~  +  tan2°eEl  "  a  tan  °eEl  km 

(33) 

where  a  is  from  (16),  The  +  choice  in  (33)  is  m<?de  such  that  (32)  yields  its 

«> 

smallest  positive  value.  If  and/or  0.^  are  not  specified,  thev  mav  be  esti- 

mated  [26,  Sec.  2.4]  using  Ah  [km]  which  is  a  model  input  parameter  [20,  p.  101]  and 
d( sl  from  (27) .  That  is: 

hg  =  larger  of  (hel  or  0.005)  km 

(34) 

dh  =  dLsl  exp  (~°* 07  /  Ah/he)  km 

(35) 

0.1 

dLsl 

if 

dh  < 

0.1 

dLsl 

dLEl  = 

3 

dLsl 

if 

dh  > 

3 

dLsl 

km 

(36) 

dh 

otnerwise 
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I 


Facility 
Horizontal 


Note:  Effective  earth  (straight  ray) 
geometry  Is  illustrated  for  the  facility 
using  a  dotted  ray.  The  solid  ray 
Illustrates  that  the  horizon  ray 
obtained  by  ray  tracing  from  the 
facility  hori2on  to  the  aircraft 
yields  a  smaller  distance 
than  would  be  obtained  with 
effective  earth  geometry. 


Not  drawn 
to  scale 


FIGURE  4.  FACILITY  RADIO  HORIZON  GEOMETRY. 
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lesser  of 


or 


rad 


(37; 


j  0.2094  (12°)  J 

where  d  is  from  (27). 

However,  some  of  this  flexibility  must  be  sacrificed  when  the  facility  is  high 
(airborne)  since  the  accurate  specification  of  more  than  one  horizon  parameter  re¬ 
quires  prior  knowledge  of  ray-tracing  results. 

The  IF-77  version  was  constructed  to  retain  all  of  the  IF-73  facility  horizon 
specification  flexibility  for  low-facility  antennas  and,  yet,  allows  ray  tracing  to 
be  used  for  high-facility  antennas.  This  method  is  described  in  the  following  steps 

1.  Determine  horizon  parameters  as  they  were  determined  in  IF-73;  but,  con¬ 
sider  the  results  as  initial  values  that  may  be  changed  if  the  facility  antenna  is 
too  high. 

2.  Values  of  from  (21)  and  hgl  from  (24)  are  used  to  test  the  initial 
horizon  parameters.  The  initial  parameter  values  are  replaced  by  ones  appropriate 
for  a  smooth  earth  if  the  test  conditions  are  met;  i.e.,  smooth-earth  values  are 
used  if; 


hel  >  3  km  and 


0 

etl 

>  0 

and 

>  hLEl 

or 

GeEl 

<  0 

and 

A 

— j 

1 — 
m 

3.  Step  3  is  not  used  if  smooth  earth  parameters  were  selected  in  step  2. 

If  Ahgj  from  (25)  is  zero,  the  initial  horizon  parameter  values  from  step  1  are  used 
Otherwise,  r«y  tracing  is  used  to  determine  values  for  0,  and  d,  •, .  In  either  case: 


eeEl  if  flhel  *  0 


otherwise  use  ray  tracing 


dLEl  if  Ahel  =  0 


otherwise  use  ray  tracing 


The  ray  tracing  referred  to  in  (40)  and  (41)  is  started  at  the  horizon  elevation, 
h[_rr  w1t^  a  take~°ff  angle  of  -6^  and  continues  until  the  facility  antenna  height 
hrj  is  reached.  Then,  the  great-circle  distance  traversed  by  the  ray  is  taken  as 
djj ;  and  the  negative  of  the  ray  arrival  angle  is  taken  as  0  , .  The  take-off  angle 
used  is  calculated  from: 


°L  “  °eEl  + 


(dLEl/a) 


Figure  5  provides  a  summary  of  the  logic  used  for  facility  horizon  determination. 

The  distance  dLR2  shown  in  Figure  4  is  taken  as  the  distance  under  a  ray  traced 
from  the  facility  horizon  with  a  take-ott  angle  of  0^  from  (42)  to  the  aircraft 
altitude  of  h^.  This  distance  is  then  used  with  dLs1  2  from  (27)  to  calculate  the 
maximum  1 ine-of- sight  distance  dM1  which  is  also  shown  in  Figure  4;  i.e.: 


dLsl  +  dLs2  ^or  sm00t^  eerth;  i.e..  Ah  =  0 


(-'  (%-! 


el  Jcos0el 
+  h”o) 


dLl  +  dLR2 


otherwi se 


°el  i  if  Ahe2  =  0  f  km 


3.3  Aircraft  (Or  Higher  Antenna)  Horizon 

Aircraft  horizon  parameters  are  determined  using  either  (a)  case  1,  where  the 
facility  horizon  obstacle  is  assumed  to. provide  the  aircraft  radio  horizon,  or  (b) 
case  2,  where  the  effective  reflection  surface  is  assumed  to  provide  the  aircraft 


1C 


FIGURE  5.  LOGIC  FOR  FACILITY  HORIZOft  DETERMINATION. 
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with  a  smooth-earth  radio  horizon.  The  great-circle  horizon  distance  for  the 
aircraft,  is  calculated  using  the  parameters  shown  in  Figure  6  along  with  the 

great-circle  distance,  d,  between  the  facility  and  the  aircraft;  i.e.: 


'sL  =  v/  2a  hL r  1 

km 

!LM  =  dLl  +  dsL 

+  dLs2  km 

d-du 

if  dML  <  d  <  d 

L2  = 

dLs2 

otherwi se 

Here,  \rl  the  height  of  the  facility  horizon  obstacle  above  the  effective  reflec 
tion  surface  from  Figure  5,  and  d^  is  the  smooth-earth  horizon  distance  for  the 
obstacle  (i.e.,  a  is  from  (16),  is  from  (41),  and  dL$2  is  from  (27)).  The 
horizon  ray  elevation  angle  at  the  aircraft,  0e2,  is  measured  relative  to  the  hori¬ 
zontal  at  the  aircraft,  wi th  positive  values  assigned  to  values  above  the  horizontal 
it  is  calculated  from: 


Horizontal  at  aircraft 
Case  1,  obstacle  horizon 
Case  2,  smoothr-earth  horizon 


Not  drawn 
to  scale 


FIGURE  6.  GEOMETRY  FOR  AIRCRAFT  RADIO  HORIZON 


(47) 


\rl  if  dHL  "  d  <  dLM 


Lr2 


km 


0  otherwise 


J 


hL2  =  hLr2  +  hr 


km 


(48) 


ee2  =  Tan 


( 


-1  [  hLr2  ‘  he2  dL2 


L2 


et 


) 


rad 


(49) 


where  h^  is  the  aircraft  horizon  height  above  the  reflecting  surface. 


4.  DIFFRACTION  REGION 

Calculations  based  on  diffraction  mechanisms  are  used  both  within  and  beyond 
the  radio  horizon.  Diffraction  attenuation,  A^,  is  assumed  to  vary  linearly  with 
distance  in  the  diffraction  region  when  other  parameters  (heights,  etc.)  are  fixed. 
Most  of  the  equations  given  in  this  section  are  related  to  the  determination  of  two 
points  needed  to  define  this  diffraction  line.  Since  irregular  terrain  may  be 
involved,  rounded-earth  diffraction  is  combined  with  knife-edge  diffraction  consid¬ 
erations.  This  is  done  by  comDining  attenuation  values  obtained  via  rounded  earth 
with  those  obtained  using  knife-edge  diffraction  at  two  distances  (d^  and  d^),  and 
fitting  a  straight  attenuation  versus  distance  line  to  them.  The  paths  involved  may 
be  illustrated  using  the  points  shown  in  Figure  7  as  path  F-O-ML  for  d^  and  path 

ML 


19 


F-O-A  for  d^  where  corresponds  to  the  maximum  1 ine-of-sight  distance  and  d^  is 
the  shortest  heyond-the-horizon  distance  that  involves  both  the  facility  horizon 
obstacle  and  a  smooth-earth  horizon  for  the  aircraft. 

Rounded-earth  and  knife-edge  diffraction  calculations  are  discussed  in 
Sections  4.1  and  4.2,  respectively.  Section  4.3  deals  with  the  determination  of  the 
diffraction  attenuation,  A^. 
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M 


P 


G_ 

hpl  ,2 


I 


Attenuation  line  slope 
Mp  for  path  F-O-ML 
Mq  for  path  0-A 

Height  gain  function 

\f1  2  ^or 

GFo]  ^  f°r  P«th  0_A 


> 

j 


dB/km 


dB 


(55) 


(56) 


With  appropriate  starting  parameters  from  (50)  through  (53),  Ap,  Mp, 
determined  as  follows: 


and  G_  are 
hpl  ,2 


a  =  effec  ve  earth  radius  from  (16) 
f  =  frequency  [MHz] 


°3  = 

0.5 

(a2/f )1/3/a 

rad 

°4  = 

393 

rad 

dLp  * 

dLpl 

+  dLp2 

km 

d3 

dLP 

+  ae3 

km 

d4 

d3  + 

2a03 

km 

‘1,2  = 

dLpl, 

2/(?hepl,2) 

km 

3,4  = 

<d3,4 

'  dLp)/03,4 

km 

o  =  conductivity  (Siemens)  which  is  an  input  parameter 
[20,  p.  99,  SURFACE  TYPt  OPTIONS] 

x  =  18000  o/f 

;  =  dielectric  constant  which  is  an  input  parameter 
[20,  pp.  99,  SURFACE  TYPE  OPTIONS] 


(57) 

(58) 

(59) 

(60) 
(61) 
(62) 

(63) 


(64; 
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R 


Kdl,2,3,4  =  °-36278  <fal,2,3.4>'1/3E(£  '  +  *‘1 


,2  ,  „2,-l/4 


Kdl  2  3  4  f°r  horizontal  polarization 


'1 ,2,3,4 


or 


2  2  1/2 

Kdl  2  3  4  +  J  ;  f°r  vertical  polarization 


Kf;i  2  =  smaller  of  (K1  2  or  0.99999) 


,1/3 


l,2,3,4=416'4f'  ('•607-K]>2,3,4) 


Xl,2=  B,,2<23dLpl,2  km 


(3,4  83,4a3^4^  (d3,4  ’  dLp*  +  X1  +  X2  km 


G1 ,2,3,4  =  °-05751  X1 ,2,3,4  _  10  109  X1 ,2,3,4 
W]j2  =  0.0134  X1  2  exp  (-0.005  X]  2) 


(65) 

(66) 

(67) 

(68) 

(69) 

(70) 

(71) 

(72) 


y1>2  =  40  log  (X^)  -  117 


dB 


1,2 


r 


y-j  2  or  -117  whichever  has  the  smaller  absolute 
value  for  0  _<  Kp^  2  <  10~5 


or 


_5 

y^  2  f°r  10  1  2  ancl 


or 


■450  ^log  Kf1>2^ 


-3 


<  X, 


i,? 


20  log  (Kn  2)  -  15  +  2.5  (10'5)X^  2/KF1 


otherwise 


(73) 


>  dB  (74) 
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*,  ..  . , 


■1  lil  rr  “  i 


1  tk  »  •  i'i  A  i  .4j 


Whenever  hcp] >2  >  2  fci>2 
0 

or  waen  2  2.  5 

-6.6  -  0.013  h[  2  -  2  log  h]  2 

or  when  2  -  °‘C^ 

1.2  -  13.5  h]  2  +  15  log  h]>2  if  h]>2  <  0.3 

-6.5  -  1.67  h1>2  +  6.8  log  h]<2  if  h]j2  >  0.3 

or  when  0.01  <  2  1  0-05 

T  -  25  (T  -  B)  (0.05  -  K^) 

«Fi1.2  “  i  where 

T  =  -13.9  +  24.1  hlj2  +  3.1  log  h12  if  h]j2  <  0.25 

T  =  -5.9  -  1.9  h1  ?  +  6.6  log  h1>2  if  h~1>2  >  0.25 

B  =  1.2  -  13.5  h1>2  +  15  log  h1>2  if  ^  ^  <  0.3 

B  =  -6.5-  1.67  h|>2  +  6.8  log  h1>2  ifh]>2>0.3 

or  when  K-j  2  > 

T  -  20  (T  -  B)  (0.1  -  K1>2) 
where 

T  =  -13  if  h-,  2  <  0.1 

T  =  -  4.7  -  2.5  h]>2  +  7.6  log  h]  2  ifh1>2>0.1 

B  =  -13,9  +  24.1  h1  2  +  3.1  log  h1>2  if  h]>2  <  0.25 

,  r  *  ...  t-  ’  >  0.25 

B  =  -  5.9  -  1.9  h-j  2  +  6.6  log  h-,  ^  if  h^2  - 
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Note  that  many  of  the  variables  used  in  (60)  through  (84)  have  values  that  are 
dependent  upon  path  parameters  from  (50)  through  (56),  but  are  not  identified  with  a 
p  subscript. 

The  formulation  provided  in  (83)  for  G^  2  is  based  on  a  curve  fit  approxima¬ 
tion  to  the  residual  height  gain  function  curves  [34,  Fig.  7.2]  developed  by 
A.  G.  Longley  of  NTIA  (unpublished  paper,  “Calculation  of  Transmission  Loss  for  Fre¬ 
quencies  from  200  MHz  to  10  GHz").  While  this  approximation  may  yield  incorrect 

values  for  K,  2  >  0.1,  h^  2  <  0.01  km  or  h^  2  >  1 00  km,  the  application  of  the 
height-gain  function  'in  many  such  cases  should  be  tempered  anyway.  Although  special 
consideration  for  the  K-j  2  >  0. 1  and  2  <  0-01  km  cases  are  not  included  in  (83), 

the  initial  test  of  hepl  2  against  2fcl  2  will  result  in  Ghl,  2  =  0  for  almost  all 

cases  where  h-|  2  >  100  km.  Furthermore,  the  abrupt  G^  2  change  associated  with 
this  test  is  tempered  by  the  "blending  function"  of  (82)  in  (84)  [24,  eq.  6]. 

Rounded-earth  diffraction  attenuation  for  the  path  F-0-A,  Arp,  is  given  by: 

Ar[-  =  Ap  +  Mpd  dB  (85) 


and  the  rounded-earth  diffraction  attenuation  for  path  0-A  at  distance  dp^  +  dpg2, 


Ar0.  is: 


Ar0  '  A0  +  H0*dL01  +  dL0?' 


dB 


(86) 


4.2  Knife-Edge  Dif Traction 

Knife-edge  diffraction  attenuation  by  an  isolated  obstacle  with  ground 
reflections  [24,  Sec.  3.5;  34,  Sec.  7.2]  is  calculated  for  the  F-0-M  and  F-0-A 
paths  [34;  Sec's.  7.2,  III. 3]  illustrated  in  Figure  7;  i.e.,  AKML  These  calcu¬ 
lations  utilize  G^-Fi  2  and  G^  from  (56)  and  (84),  dp(^  2  from  (52,53),  0^  and 
d^  from  Figure  5,  0es2  from  (30),  a  from  (16),  ArQ  from  (86),  dMp  from  (43),  and 
frequency  f[MHz].  That  is: 


AKML  =  6  "  GhFl  “  GhF2 


dB 


(87) 


dLsA  =  dL01  +  dL  02 


km 


(88) 
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0v  ■  eel  +  °es2  +  <dU  *  W'8  rad 


(89) 


v  =  5.1658  sin  (0.50^) 


Cv  =  / q  cos(-|  )  dt,  Sv  -  Jq  sin  )  dt 


where  these  are  Fresnel  integrals  [34,  p.  III-18]. 


fv=0.5  v[1  -  (Cy  +  Sy)]2  +  (Cw  -  Sw) 


v  V‘ 


AKA  ArO  "  ^hFl  ‘  6h01  "  20  log  fv  dB 


4.3  Diffraction  Attenuation, 

Diffraction  attenuation,  Ad,  is  calculated  using  the  rounded-earth  diffraction 
(Sec.  4.1)  and  knife-edge  diffraction  (Sec.  4.2)  parameters  just  discussed;  i.e. 

/ 

1  when  d^L  _>  dps  (rounded-earth  only) 

0  when  d^  £  0.9  dps  (knife-edge  only) 


0.5(,+cos(%^)) 


otherwise 


(combination  of  both) 


where  the  maximum  1 ine-of-sight  distance,  d^ ,  is  from  (43),  and  the  total  smooth- 
earth  horizon  distance,  d.  ,  is  from  (28): 


f  ArML  lf  W  ■'  °*! 

/  A  •  /-  la  i 


akml  if  w  "  °-001  dB 

(1  -  W)  Akml  +  WArML  otherwise 


where  ArML  =  ArF  from  (85)  with  d  =  dM  and  l\  is  from  (87) 


dA  =  d.  ,  +  d.  .  kin 

A  LI  LsA 


"  VV  V*  V  . V* . 


is  from  Figure  5,  and 


is  the  faci 1 i ty-to-aircraf t  distance  of  figure  7,  d^ 
^LsA  ^rom 


{ArA  if  W  0.999 
A^  if  W  <  0.001 
(1  -  W)  A^  +  WAp  otherwise 


dB 


where  AfA  =  A,,p  from  (85)  with  d  =  d^  and  A^A  is  from  (93). 


Md  "  <AML  -  V^ML  "  dA>  dB/km 


■  v  ■ 

Md  dML 

dB 

=  fld0  + 

Mdd 

dB 

(97) 


(98) 

(99) 
(100) 


where  d  fkm]  is  the  great-circle  path  distance  and  A^  is  the  diffraction  attenuation 
applicable  to  it.  The  use  of  Ad  in  the  determination  of  terrain  attenuation,  Ap 
is  discussed  in  Section  7. 


5.  LINE-0F-SIGIIT  REGION 

Calculations  based  on  a  two-ray  model  are  used  in  estimating  an  attenuation, 
AL0S  ^ec.  5.4),  and  short-term  fading  statistics,  (Sec.  10.2),  for  paths  with 
distances  less  than  the  maximum  1 ine-of-sight  (LOS)  range,  d^,  from  (43).  This 
model  involves  the  phasor  summation  of  the  earth  rc  ^cted  and  direct  rays.  The 
geometry  and  effective  reflection  coefficient  associated  with  it  are  discussed  in 
Sections  5.1  and  5.2,  respectively.  At  distances  less  than  d^p  but  greater  than 
the  LOS  transition  distance,  dQ  (Sec.  5.3),  diffraction  makes  a  contribution  to 

alos* 


5.1  Two-Ray  Path  Length  Difference  Geometry 
Geometry  for  the  two-ray  model  used  in  the  I  OS  region  is  shown  in  Figure  8. 
Path  length  difference,  Ar,  is  the  extent  by  which  the  length  of  the  reflected  ray 
path,  r]  +  r2  =  ri2>  exceecls  that  of  the  direct  ray  rQ;  i.e.: 
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km 


001) 


The  geometrical  parameters  encountered  in  the  calculation  of  A r  are  used  in  other 
parts  of  the  model  (e.g.,  antenna  gain  factors.  Sec.  5.2.4),  and  Ar  is  used  in  the 
two-ray  phasor  summation  that  results  in  the  LOS  lobing  structure  (Sec.  5.4). 

The  Ar  formulation  provided  here  uses  aQ  from  (15),  a  from  (16),  Ahgl  2  from 

(25),  h  ,  9  from  (20),  the  height  of  the  counterpoise  above  ground  h  ,  and  the 
r  I , c  eg 

height  of  the  facility  antenna  above  its  counterpoise,  h^c«  Both  hc^  and  h^c  are 
model  input  parameters  and  are  used  only  if  the  facility  counterpoise  option  is 
desired.  This  formulation  provides  the  Ar  and  path  length,  ci,  for  a  particular 
grazing  angle,  ^ ;  and  it  is  repeated  for  various  values  of  ip  until  sufficient 
coverage  ’ ~  obtained  for  both  Ar  and  d.  If  a  counterpoise  is  present,  a  set  of  cal¬ 
culations  will  also  be  required  with  the  counterpoise  parameters  selected.  The 
formulation  may  be  summarized  as  follows: 


(102) 

(103) 

(104) 

(105) 

(106a) 

(106b) 

(107) 

(108) 

(109) 
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u  ^  M,  -  -  |«.  ™-  mj>  m  •  my« 


Hi,2 


^  tan  V  f°r  <!'  <  1-56  rad  *1 


2  otherwise 


) 


km 


(110) 


a 


Tan_1[(H^  -  Hp/(D]  +  D2)]for  ip  <  1.56  rad  1 


ip  otherwise 


km 


(111) 


ro 


12 


(D^  +  02)/cos  a  for  i p  <  1.56  rad 


otherwise 


(D1  +  D2)/cos  ^  for  <  1.56  rad 


H1  +  H2  otherwise 


km 


km 


Ar  =  4  HJH£/(r0  +  r]2)  km 


Ary  =  Ar  with  earth  parameters  in  (106)  km 


Ar  =  Ar  with  counterpoise  parameters  in  (106)  km 

V 


ehl  =  a  -  6^  rad 


0h2  -  -(a  ♦  e2) 


Vl  ,2  '  +  °1,2* 


rad 


rad 


0  =  0,  +  0o  rad 

o  1  2 


d  =  a  0  km 
a  o 
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(112) 


(113) 

(114a) 

(114b) 

(114c) 

(115a) 

(115b) 

(116) 

(117) 

(118) 


An  adjusted  effective  earth  radius,  aa,  and  adjusted  effective  antenna  heights, 
2*  that  vary  with  i|>  are  used  in  this  formulation  since  the  values  of  a  from  (16), 
and  hrj  ^  from  (20)  are  not  appropriate  for  large-ray,  take-off  angles  when  cos  ip  is 
not  near  1  [3,  eq.  3.23],  Since  rQ  and  r^  can  be  very  large  and  nearly  equal,  Ar 
is  calculated  via  (114)  Instead  of  (101).  Except  when  the  counterpoise  case  is 
specifically  mentioned,  future  references  to  variables  calculated  in  (106)  through 
(118)  involve  the  earth  reflection  case  (e.g.,  Ar  =  Ar  ) . 


5.2  Effective  Reflection  Coefficient' 

A  counterpoise  for  the  facility  antenna  (e.g.,  as  with  very  high  frequency 
omni-range  or  VOR  antenna)  is  an  option  of  the  1F-77  model  [20,  p.  89].  Hence, 
reflection  coefficients  for  both  earth  and  counterpoise  reflections  are  computed. 
Magnitudes  for  these  coefficients  are  and  R^c,  respectively.  The  factors  in¬ 
volved  in  these  coefficients  include  the  divergence,  Dv,  and  ray-length  factors,  F 
(Sec.  5.2.1);  the  surface-roughness  factor,  Fa^  (Sec.  5.2.2);  counterpoise  factors, 

f  .  (Sec.  5.2.3);  antenna  gain  pattern  factors,  gD„  „  (Sec.  5.2.4);  and  plane-earth 
g  *C  KQ  *c 

reflection  coefficient  magnitudes,  Rcg  (Sec.  5.2.5);  i . e . : 

RTg  =  WohVRgRg  °'9a) 


0  when  there  is  no  counterpoise 


fcgRcRc  otherwise 


J 


(119b) 


Since  the  counterpoise  is  taken  tc  be  near  the  facility  antenna  and  flat,  DyFrFah 
is  taken  as  unity  in  (119b).  It  there  is  no  counterpoise  f  =  1 ,  Ryc  =  0,  and  the 
calculation  of  f  ,  gD  .  and  R_  is  not  done.  When  isotropic  antennas  are  used, 

C  KC  C 

gRg,c  =  L 


5.2.1  Divergence  and  Ray  Length  Factors 
Reflection  from  the  curved  earth  surface  is  less  efficient  than  reflection  from 
a  flat  earth  would  be  [15,  Sec.  3.2;  35).  This  reduction  is  taken  into  account  by  a 
divergence  factor,  D  ,  which  multiplies  the  plane-earth  reflection  coefficients  as 
in  (119a). 
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When  both  antennas  are  high  and  close  (i.e.,  two  close  aircraft),  the  relative 
magnitude  of  the  reflected  ray  will  be  reduced  because  the  reflected  ray  length  is 
much  longer  than  the  direct  ray  length  (i.e.,  larger  free  space  loss  for  reflected 
ray)  [15,  Sec.  3.3].  The  ray  length  factor,  F  ,  in  (119a)  is  used  to  account  for 
this. 

The  formulation  for  Dy  and  F  may  be  summarized  as  follows: 


rT,2 


J  H]  2  if  <p  s  90° 

2/cos  ?  otherwise 


knr 


(120) 


where  H^2  are  from  (106);  the  grazing  angle,  (Fig.  8),  is  a  starting  parameter 
for  the  Ar  formulation  of  Section  5.1;  and  D1  2  are  from  (109) 


Rr 


rlr2 

r12 


km 


(121) 


where  r]z  are  from  (120)  and  r]2  is  from  (113) 


Dv  = 


1  + 


2Rr(l  +  sin  ip) 
a  sin-? 

a 


ft)'] 


-1/2 


(122) 


where  afl  is  from  (104),  and 


Fr  =  ro/r12 


where  rQ  is  from  (112). 


(123) 


5.2.2  Surface  Roughness  Factors  ; 

Surface  roughness  factors  for  specular,  F  and  diffuse,  F^  reflections 
are  calculated  as  follows: 

•\hd  =  Ahm[l  -  0.8  exp  (-0. 02d ) ]  in  (124) 
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where  Ah  is  the  terrain  parameter  used  to  characterize  irregular  terrain  and  is  an 
IF-77  model  input  parameter  [20.  P-  101];  and  d  is  great-circle  path  distance  from 

(118). 


0.25  H1/c 
0.39  Ah, 


for  water 


for  Ah^  £.  4  m 


0. 78  Ahj  exp  ^-0. 5 


Ahd'/4 J otherwise 


(125) 


where  is  a  value  for  significant  wave  height  that  is  selected  based  on  the  sea 

state  when  a  water  reflecting  surface  option  is  chosen  [20,  p.  99]. 


6  =  o.  sin  (i|i)/A  m 
h  m 


(126) 


where  the  grazing  angle,  i p  (Fig.  8),  is  a  starting  parameter  for  the  Ar  formulation 
of  Section  5.1;  and  the  wavelength  in  meters,  Am,  is  from  (10a). 


Fq|1  =  exp  (-2ir6) 


(127) 


0.01  +  94662  if  6<0. 00325 
6.156  if  0. 00325<6<0. 0739 

0.45  +  v/0. 000843  -  (6  -  0.1026p  if  0.0739'6<0. 1237 

0.601  -  1.066  if  0. 1 237<6<0. 3 

0.01  +  0.875  exp(-3.886)  otherwise 


The  formulation  for  Fdoh  is  based  on  curves  fit  to  data  [5,  Fig.  4],  and  is 
used  along  with  the  plane-earth  reflection  coefficient  magnitude,  Rg  (Sec.  5.2.5); 
ray  length  factor,  Fr,  from  (123);  counterpoise  factor,  f  (Sec.  5.2.3);  and  the 
antenna  gain  pattern  factor,  gR^  (Sec.  5.2.4);  to  formulate  the  diffuse  reflection 
coefficient;  i.e.: 


Fr  ^doti  fg  ^Rg  % 


i 


■  ■;* 


■  v’-k  *  ; '  w  '-  .  >  ,■>  v  *vr*  ’ 


5.2.3  Counterpoise  Factors 

Counterpoise  factor,  f  ,  in  (119a)  is  used  to  provide  some  reduction  in  the 
reflection  from  the  earth's  surface  when  this  reflecting  surface  is  shadowed  by  the 
counterpoise.  Counterpoise  factor,  f  ,  in  (119b)  is  used  to  provide  some  reduction 
in  the  reflection  from  the  counterpoise  because  of  the  limited  area  of  the  counter¬ 
poise.  When  there  is  no  counterpoise,  f  =  i  and  f  =  0.  These  factors  are  deter- 

9  c 

mined  with  the  geometry  shown  in  Figures  9  and  10  by  using  knife-edge  diffraction 
considerations  as  follows: 


9ce  "  Tan  (2  hfc/dc) 


(130) 


where  h^c  is  the  height  of  the  facility  antenna  above  its  counterpoise  from  (22), 
and  dc  is  the  counterpoise  diameter  which  is  an  IF-77  model  input  parameter 

[20,  p.  88]. 


rc  =  0.5  dc/cos  0ce 


(131) 

(132) 


where  0^  is  from  (116): 


»,  =  rsrji 


where  A  is  from  (10b) 


0.  =  0  0.  , 

kc  ce  hi 


(133) 


(134) 


where  0^  is  calculated  using  counterpoise  parameters  in  (106)  through  (118) 


v9  *  iyv  sin  (V2) 


-  for  l°rl r  °ce 
+  otherwise 


(135a) 


vc  =  +  Vv  sin  <"kc/2) 


for  “hi  '"ce 


+  otherwise 


(135b) 


Earth 


FIGURE  9.  GEOMETRY  FOR  DETERMINATION  OF  EARTH  REFLECTION  DIFFRACTION 

PARAMETER,  v  ,  ASSOCIATED  WITH  COUNTERPOISE  SHADOWING. 
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Image 


FIGURE  10.  GEOMETRY  FOR  DETERMINATION  ,'F  COUNTERPOISE  REFLECTION 
DIFFRACTION  PARAMETER,  v  ,  ASSOCIATED  WITH  THE  LIMITED  REFLECTING 
SURFACE  OF  THE  COUNTERPOISE. 


■,r-: 


v" 


Cq  c  =  /  9,C  cos 
9,c  q 


M 


dt  ,  S  =  f 

g,c  J 


;-c 


where  these  are  Fresnel  integrals  [34,  p.  III-18], 

fg,c  =  °-5V[l  -  (Cg,c  +  ,Sg,c^2  +  (C9>c  -  sg,c): 


kg,c 


=  Tan 


-1 


v  ”  ug,c  "  Sg,c  / 


rad 


(136) 


(137) 


(138) 


The  angles  <j>.  are  phase  shifts  that  will  be  used  later  in  Section  5.4. 
Kg ,  c 


5.2.4  Antenna  Pattern  Gain  Factors 

The  ar.tenna  gain  factors  gQ  p  and  g^  v  are  used  to  allow  for  situations  where 
the  antenna  gains  effective  for  the  direct  ray  path  differ  from  those  for  the 
reflected  ray  path.  Figure  11  illustrates  the  two-ray  path  and  indicates  the  gains 


Note:  This  sketch  is  drawn  with  flat  earth,  straight  rays  and  an 
exaggerated  scale  so  that  the  geometry  shown  is  over  simplified. 


FIGURE  11.  SKETCH  ILLUSTRATING  ANTENNA  GAIN  NOTATION 
AND  CORRESPONDENCE  BETWEEN  RAY  TAKE-OFF  ANGLES  AND  GAINS. 


$ 


g/.IV'.i'.iv!;. 
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involved.  These  gains  are  the  relative  voltage  antenna  gains  (volts/volt  or  V/V). 
They  are  measured  relative  to  the  main  beam  of  their  respective  terminal  antennas; 
i.e.,  for  main  beam  conditions  g^  R=gR^  £=1  V/V.  This  convention  is  consistent  with 
usage  in  IF-73  [14,  p.  39].  However,  it  is  NOT  CONSISTENT  with  usage  in  the  Multi- 
path  Handbook  [16,  Sec.  CI-D.3]  where  identical  symbols  are  used;  but,  the  gains  are 
measured  relative  to  an  isotropic  antenna. 

In  general,  these  gains  are  complex  quantities;  but,  IF-77  includes  provisions 
for  scalar  gains  only;  i.e.,  these  gains  are  >0  and  <] .  In  many  practical  applica¬ 
tions,  the  direct  and  reflected  rays  will  leave  (or  arrive)  at  elevation  angles  where 
the  relative  phase  is  either  expected  to  be  near  zero  or  is  unknown  so  that  the 
complex  nature  of  these  gains  is  largely  academic.  They  are  called  voltage  gains 
since  they  are  a  voltage  ratio  that  could  be  considered  dimensionless  (volt/volt), 
but  are  different  from  gains  expressed  as  power  ratios  (watt/watt)  that  could  also 
be  considered  dimensionless.  Decibel  gains  above  main  beam  values  are  related  to 
these  gains  by  formulas  such  as: 

GR1 ,2^dB^  =  20  log  gRl,2  (139) 

where  GR1  2  <  0  because  0  <  gR1  2  <  1 
and 

/Gfj-i  o/20) 

gRl,2[V/V]  =  10  ’  (140) 

The  formulations  for  gD  R  are: 
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where  isotropic  implies  that,  for  the  radiation  angles  of  interest,  gD1  =  g.,..  and 
gR2  =  9q2’  In  Problems  "involving  elliptical  polarization,  horizontally  polarized 
(ghDl,2  and  ghRl,2^  and  vertical1y  polarized  (gyD1  2  and  gyR^  2)  components  are  used. 
Linear  polarization  is  considered  to  be  either  vertical  or  horizontal  with  the  polar¬ 
ization  associated  with  gQ  R  selected  accordingly.  Defining  gR  as  1  for  elliptical 
polarization  is  done  to  allow  the  antenna  gains  to  be  included  in  the  reflection 
coefficient  formulation  of  IF-77  .n  a  simple  way  for  horizontal  or  vertical  polari¬ 
zation.  Circular  polarization  is  a  special  case  of  elliptical  polarization;  i.e., 


ghDl,2  9vD1 ,2 


The  gain  factor  gRy  is  similar  to  gR  except  that  gRv  involves  gains  gyR1  2;  i.e. 


9rv[V7V] 


1  for  isotropic  antennas 
gvRlgvR2 


T 

otherwise  j 


In  a  similar  manner. 


tRh[V/V]  ■ 


1  for  isotropic  antennas 


^ghRlghR2  otherwise 


M43) 


(144) 


where  gRh  is  for  horizontal  polarization.  These  factors  will  be  used  in  the  formu¬ 
lation  of  complex  plane-earth  reflection  coefficients  for  elliptical  polarization 
that  is  given  in  the  next  section. 

Several  facility  antenna  patterns  from  which  gain  factors  can  be  determined 
are  included  in  the  computer  programs  that  utilize  IF-77  [20,  p.  85].  However, 
data  for  other  facility  antenna  patterns  or  aircraft  antenna  patterns  can  be  used. 
These  programs  also  include  an  option  to  tilt  the  main  beam  of  either  or  both 
antenna(s)  relative  to  the  horizontal,  or  have  either  or  both  antenna(s)  track  the 
other  with  its  main  beam  [20,  p.  89].  The  patterns  involved  here  are  vertical  plane 
antenna  patterns.  Gain  variations  with  aznnuth  can  be  accommodated  by  adjusting 
Gj  R  in  (4),  and  R  for  (4)  are  obtained  from: 
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20  log  if  facility  is  transmitting 


(145a) 


20  log  gD2 
20  log  gD1 
(^20  log  gD2 


otherwise 

if  facility  is  receiving 
otherwise 


(145b) 


The  elevation  angles  to  be  used  in  the  determination  of  the  various  gain 
factors  are  listed  below: 


FACTOR 

gDl 

0H1 

gD2 

0H2 

gRl*gvRl ,ghR1 

egi 

gR2’gvR2’ghR2 

eg2 

ELEVATION  ANGLE 

for  LOS,  otherwise  6^  (Fig.  5) 
for  LOS,  otherwise  0g2  from  (49) 


These  angles  are  calculated  as  follows: 


(146) 


where  Qesl  2  is  from  (30),  9^  2  is  from  (29),  afl  is  from  (104),  aQ  is  from  (15), 
and  a  is  from  (16) , 


HI, 2  9hl ,2  +  eLl,2 

rad 

gl  ,2  =  0rl,2  +  0L1,2 

rad 

(147) 

(143) 


where  0^  2  is  from  (115),  and  0^  2  is  from  (116).  The  effect  of  0^  2  is  to  force 
0^  2  and  0  ^  2  to  have  the  values  obtained  via  ray  tracing  at  the  smooth-earth 
horizon,  and  prorate  values  obtained  elsewhere. 


When  a  counterpoise  is  present,  it  will  have  a  set  of  gain  factors  associated 
with  it  where  the  earth  or  counterpoise  set  values  are  determined  by  the  values  used 
for  Hj  2  1n  (106)  to  calculate  0^  ^  anc>  er]  2*  fr°r  examP^e: 


_  calculated  with  parameters  appropriate  for  a 
gRg  ~  9P.  ground  reflection  (149a) 


_  calculated  with  parameters  appropriate  for  a 
3Rc  "  gR  counterpoise  reflection  (149b) 


5.2.5  Plane-Earth  Reflection  Coefficients 
Values  for  the  plane-earth  reflection  coefficient,  R  exp(-j<J>),  are  dependent 
upon  the  dielectric  constant,  e,  and  conductivity,  a,  of  the  surface  involved.  For 
water: 


f2T(e  -  e  ) 

o[mho/mj  -  - ^§^5 -  +  ai  (151) 

where  is  the  static  dielectric  constant,  =  4.9  is  the  dielectric  constant 
representing  the  sum  of  electronic  and  atomic  polarizations,  f[MHz]  is  frequency, 
T[ps]  is  relaxation  time,  and  a..[mho/m]  is  the  ionic  conductivity.  The  values  for 
e$jT,  and  for  water  [15,  p.  26]  were  obtained  using  Saxton  and  Lane  [36].  When 
(150,151)  are  not  used,  appropriate  e  and  o  values  are  taken  from  the  Applications 
Guide  [20,  p.  89].  The  formulation  for  R  exp  (-  j$)  may  be  summarized  as  follows: 

ec  =  e  -  j  60  Xm  0 

where  A„  is  from  (10a). 
m 


(152) 


where  4>  (Fig.  8)  is  a  starting  parameter  for  the  Ar  formulation  of  Section  5.1. 

e  sin(^)  -  V 

Rv  exp:-J(n-cv)]  =  sin(l)l)  +  Y‘  9r  054a) 

where  gR  is  from  (142). 

sin(«|»)  -  Y 

Rh  e*Pt-j("-ch)]  >  sTn^)  +  ^  9„  (154b) 

Re  exp[-j(TT-ce)]  =  0.5  (gRhRh  exp[-j(7T  -  ch)]  +  gRyRy  exp  [-j(ir  -  cy]}  (154c) 


where  gRv>h  are  from  (143,144),  and 


R  exp(-j<j>) 


< 


< 


Ry  exp  [-j(tt  -  c  )]  for  vertical  polarization 
(electric  field  in  plane  of  incidence) 

Rh  exp  C-j ( 7r  -  ch)]  for  horizontal  polarization 
(electric  field  normal  to  plane  of  incidence) 
Rg  exp  [ - j ( tt  -  cg)  for  elliptical  polarization 


> 


(155) 


The  part  of  this  formulation  for  elliptical  polarization  is  valid  only  when  the 
transmitting  and  receiving  antenna  have  the  same  sense  and  circular  polarization 
is  a  special  case  of  elliptical  polarization. 

When  earth  (or  ground)  reflection  parameters  (i.e.,  gain  factors  and  surface 
constants)  are  used  in  (150)  through  (155),  the  resulting  reflection  coefficient 
is  taken  as  Rg  exp  ( -  j  ) ;  and,  when  counterpoise  parameters  are  involved,  the 
resulting  coefficient  is  R  exp(-j4>  ).  That  is: 

«  V 


Rgexp(-j<}>g)  =  R  exp(-j<j>) 


when  gain  factors  and  surface  constants  . 

appropriate  for  a  ground  reflection  are  used  '  a' 


when  gain  factors  and  surface  constants 

Rrexp(-j^  )  =  R  exp(-j(fi)  appropriate  for  a  counterpoise  reflection  (156b) 

are  used 
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where  the  geometry  used  in  calculating  gain  factors  (Sec.  5.2.4)  is  dependent  on 
the  use  of  H-|  ^  values  in  (106)  that  are  appropriate  for  either  ground  or  counter¬ 
poise  reflection. 

The  total  phase  lag  of  the  reflected  ray  relative  to  the  direct  ray  for  ground, 

t 

or  counterpoise,  4>jc,  is  given  by: 


pTg,c 


( 2-rrAr  JX)  +  <t>  +  „  +  nv*  J 2 


g,c 


rg,c 


kg,c 


g»c 


rad 


(157) 


where  Ar  is  from  (114),  X  is  from  (10b),  <j».  „  is  from  (138),  and  v„  „  is  from 

y » ^  *g « c  g,c 

(135).  If  there  is  no  counterpoise,  the  last  two  terms  of  (157)  may  be  neglected 
since  they  are  the  phase  lag  introduced  by  knife-edge  diffraction  over  the  counter¬ 
poise. 


5.3  Line-of-Sight  Transition  Distance,  dQ 
The  largest  distance  in  the  1 ine-of-sight  region  at  which  diffraction  is  con¬ 
sidered  negligible  is  dQ.  In  the  IF-73  model,  it  was  estimated  using  the  distance 
at  which  the  attenuation  associated  with  a  modified  diffraction  line  is  zero 
[14,  p.  66].  The  IF-73  dQ  is  called  d^  and  in  the  IF-77  model  is  calculated 
as  follows: 

0h5  =  2  Sin  1  [(sTTtfe)  ydML/tfdLl^dML  '  dLl^  ]  rad  (158* 
where  dML  is  from  (43),  f [MHz]  is  frequency,  and  dL1  is  from  Figure  5: 

05  =  eh5  -  Tan'1{[(h1  -  hL1)/du]  -  du/2a}  (159) 

where  hL1  is  from  Figure  5,  h^  is  from  (21),  and  a  is  from  16. 

hem2  =  h2  "  Ahe2  km  060) 

where  is  the  aircraft  altitude  above  msl  (Fig.  2),  and  Ah^  is  from  (25). 
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dL5  =  "a  tan  ^5)  tan  Q5)2  +  2a(h£ 


where  a  is  from  (16),  and  from  Figure  5. 


d5  “  dL5  +  dL1 


_  T,.-l  i  hL1  '  hem2  dL5  \ 

ee5  - Ta"  V — —  -  -a- ) 


\  =  9el  +  ee5  +  <Va)  ra<i 


=  5.1658  sin  (0.5  0g)  / fd^-j d^^/dg 


Cc  =  /  cos  )dt,  Sr  =  jJ*  sin  (5^-  )dt 


where  these  are  Fr?j‘  -*1  Integrals  [34,  p.  II 1-18] . 


5  y  [1  -  (c5  +  s5)]2  +  (c5  -  s5) 


Ak5  =  -20  log  f5 


Ar5  -  Ap  +  MFd6 


where  Ap  is  from  (54,  78)  and  Mp  is  from  (55,77), 


Ar5  if  W  >  0  099 
Ag  =  AK5  if  W  <  .j| 


(1  -  W)  A,,.-  +  WA  n  otherwise 


where  W  is  from  (94),  and 


dd  =  dML  ‘  AML^dML  ‘  d5^AML  ‘  V  km  ^ 

where  A^L  is  from  (95)  and  d^L  from  (43). 

Values  estimated  for  dQ  in  IF-73  have  been  found  to  be  too  small  when  low 
antennas  are  used  for  both  antennas.  To  correct  this  difficulty,  dQ  estimates  in 
I F“ 77  are  made  using: 


do  ' 


d^  when  d^  >  dd 


dX/6  when  dX/6  >  dLl  and  dd 


km 


(172) 


V  dd 


otherwise 


where  dL1  is  the  horizon  distance  for  the  lower  terminal  (Fig.  5);  d^6  is  the  dis¬ 
tance  at  which  the  path  length  difference,  Ar, in  a  two-ray  1 ine-of-sight  formulation 
is  equal  to  A/6  (A  is  wave  length);  and  dd  is  the  dQ  of  IF-73  [14,  p.  66].  The  dis¬ 


tance  dy6  is  the  largest  distance  at  which  a  free-space  value  is  obtained  in  a 


two-ray  model  of  reflection  from  a  smooth  earth  with  a  reflection  coefficient  of  -1. 
A  value  for  d^6  can  be  determined  by  the  repetitive  use  of  the  1  ine-of-sight  formu¬ 


lation  (Sec.  5.1)  to  define  the  Ar  to  distance  relationship;  i.e.,  (102)  through  (118) 


5.4  Line-of-Sight  Attenuation,  A^ 
Line-of-sight  attenuation,  A^q<j,  is  calculated  as  follows: 


1  if  lobing  option  is  used  and  Arg 
1  if  Arg  <  0.5  A  and 

|  gp  +  RJg  exp(-j*Tg)|  <  gD 
\  0  otherwise 


<  10A 


(173) 


where  Ar  is  from  (114),  A  is  from  (10a),  gQ  is  obtained  using  earth  reflection 
geometry  as  indicated  in  Section  5.2.4,  and  RT  exp(-j4>T  )  is  from  (119,157). 


WR0  I  9D  +  FfsRTg  sxp ( - j 4>-j-g )  +  )  I  0-0001  g^  (174) 


AR0  =  10  '°9  (gD/WR0>  dB 


075) 


M 


L 


(176) 


where  AML  is  from  (95),  Aq  =  ARq  from  (175)  evaluated  for  gQ  and  WRq  values  that 
are  applicable  to  d  =  d0,  d^^is  from  (43),  and  dQ  is  from  (172),  and 


AR0  if  d  <  do 


10S 


dB 


ML  <d  "  do)  if  d0  —  d  —  dML 


(177) 


where  d  is  from  (118).  The  lobing  option  mentioned  in  (173)  allows  lobing  due  to 
ground  reflections  to  be  calculated  for  the  first  ten  lobes  inside  the  smooth-earth 
radio  horizon  [20.  p-  99].  Note  that  lobing  associated  with  reflection  from  the 
counterpoise  is  always  included,  but  that  (119b)  gives  R^c  =  0  for  (174)  when  there 
is  no  counterpoise  (Sec.  5.2). 


6.  SCATTER  REGION 

The  Rice  et  al.  [34,  Sec.  9]  method  to  calculate  attenuation  for  tropospheric 
scatter  in  IF-73  [14,  Sec.  A. 4. 4]  is  not  applicable  to  paths  that  involve  a  very 
high  antenna  such  as  a  satellite.  This  method  was  reformulated  by  Dr.  George  Hufford 
(DOC-BL,  informal  communication)  to  include  geometric  parameters  associated  with  very 
high  antennas  where  these  parameters  are  determined  using  ray  tracing  techniques. 

The  resulting  formulation  has  been  incorporated  into  IF-77  and  may  be  summarized  as 
follows: 


al  ,2 


(hLl,2  ~ 
a  +  h 


"  hel,2* 


+  2 


e,l  ,2 


-  (W2]/si"  (¥) 


rad  (178) 
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where  hLrl  2  is  from  (39,  47),  a  is  from  (16),  the  effective  reflection  surface 
elevation,  hp,  is  an  input  parameter,  and  2  is  fro,n  (41,  46). 


ds  =  d  -  du  -  dL2 


(179) 


where  the  great  circle  path  distance,  d,  may  be  taken  as  an  input  parameter. 


es  =  9al  +  ea2  +  ds/a 


■i( 


2a  +  da2j  ds  +  hL2  "  hLl 


]/• 


0„  km 
s 


(180) 


(181) 


where  h^j  2  are  from  (38,  48) 


dz2  =  ds  "  dzl 


(182) 


^  =  157  (10'6)  per 


(183) 


dN  =  A  -  a'  per  km 


Ye  =  Ns  (10”  )/dN  km 


(184) 

(185) 


where  N  is  from  (14) 


al ,2  "  7a 


k  (t*0  +  eal,2  (¥) 


+  hLrl ,2  km 


(186) 


zbl,2  =2a  dzl,2  +  0al,2  dzl,2  +  hLrl,2  km 


(187) 


Qol,2  =  Am  ”  dN  exP(-hLrl,2/Ye) 


(188) 


Qal,2  =  Am  ”  dN  exP('zal,-2^e) 


(189) 


Qbl  .2  ’  \  -  dN  exp(-zWj2/Ye) 


■al,2 


(70, 


ol  ,2  +  6Qal  ,2  ’  Qbl,2;  TS 


zK2 


+  0 


a  1 . 2  J 


]zl  ,2 


+  hLrl,2  km 


■bl,2 


^Qol ,2  +  2Qal ,2^  Z6’2+eal,2 


dzl ,2  +  hLrl ,2  km 


QA1,2  =  Am  "  dN  exP(-Zals2/Ye} 


QB1.2  =  Am  '  dN  exP(-Zbl,2/Ye) 

ZA1 ,2  =  ljQol,2  +  2QAi  ,2^  +  0al,2. 


dzl ,2  +  hLrl ,2  km 


JA1,2 


zl  ,2 


al  ,2 


rad 


e  =  eAi  +  eA2  rad 


Xa  ^ZA2  "  ZA1^9  km 


dZl  =  dzl  +  Xa 


km 


dZ2  ~  dz2  "  Xa 


km 


h''  =  (0A1ZA2  +  0A2  ZA]^/6  km 


(1 


0 


(1 


(1 


(1 

(1 

(1 


(1 


(1 


(1 


(2 

(2i 


where  h  is  the  height  of  the  scattering  volume. 


XA1,2  =  (hel,2‘hLl,2)  +  4<a  +  hel,2)(a+  hLl,2)  sin  \ 


£1,2  =  AU2  +  dZ1 ,2  km 

£,  =  +  ^2  km 

where  £  is  the  total  ray  length. 

e]  =  0.031  -  2.32  (10"3)Ns  +  5.67  (10'6)N2 

e,.  =  0.0002  N2  -  0.06  N  +  6.6 
2  S  S 

Y  ”  0.1424{1  +  e^yexp[(hv/4)^]}  per  km 

e2 

S  =  83.1 - ^ - 7 

e  1  +  0.07716  hj 

+  20  iog[(0.1424/Y)2exp(yhv)]  dB 

where  Se  is  the  scattering  efficiency  term, 
s  *  (£-|  -  £2)/£ 

where  s  is  the  modulus  of  asymmetry. 

A  =  (1  -  s2)2 
n  =  yOi/2 
Xvl  =  (1  +  s)2n 
X1(o  =  (1  -  s)2n 


where  k  is  the  wave  number  and  A  is  from  (10b). 


where  h 


where  S, 


pl,2  =  2K0hel,2 
el  2  are  from  (24). 


2  2 

ql,2  =  Xvl ,2  +  P1 ,2 


Bs  «  6  +  8s' 


+8(1+  s)X22p2/q2 
+8(1-  s)X2]P2/q2 
+  2  (1  -  s2)(l  +  2X^)0  +  2X22/q2) 


'1 


r  =  12  UL.+  +  P1  +  P2 

s  \  p 


2  /  \P]  +  p2  +  2/?> 


Sv  =  10  log 


((An2  +  B  n)q1q2  \ 

-  S--L1  +  c$  I 

plp2  f 


dB 


I  is  the  scattering  volume  term,  and 
A s  =  Se  +  SV  +  10  109  (k03/?) 


dB 


where  A  is  the  tropospheric  scatter  attenuation. 


7.  TERRAIN  ATTENUATION,  Ay 

The  IF-77  formulation  for  terrain  attenuation,  Ay,  may  be  summarized  as 
follows: 


at  * 


1 


AL0S  f0r  d  <  dML 
Ad  if  Asx  i  Adx 


A  + 
v  sx 


Asx  ~  AML 
dx  "  dHL 


for 


(d  -  dx)  otherwise!  dML<d<dx 


lesser  of  l\A  or  A,. 

d  s 

A„  once  A„  has  been  selected 
s  s 

for  a  shorter  distance 


for  dx<d 


dB  (221) 


where  ALqC  is  from  (177),  d  is  a  specified  parameter  except  in  the  LOS  region  where 
it  is  calculated  via  (118),  Ad  is  from  (100),  A$x  *  Ag  in  (220)  for  d  =  dx> 

Adx  =  Ad  ™  (10°)  for  d  ~  dx’  AML  ^rom  (95),  dML  is  ^rom  (43)»  and  As  is  from 
(220).  The  distance,  d  ,  is  the  shortest  distance  just  beyond  the  radio  horizon 

(dM^  <  d)  at  which  A$  >  20  dB  and  M$  <_  Md  where  is  the  slope  of  the  Ag  versus 

d  curve  as  determined  by  successive  Ac  calculations  (Sec.  6)  and  is  from  (98). 


FfvTL 

The  IF-77  formulation  t 
for  use  in  (3),  may  be  summar 


LOSS,  L^y 

:  ee  space  basic  transmission  loss,  Lby, 

O'fis: 


rWH  =  V  (h2  ‘  *  4<h]  +  a0^h2  +  ao)Csin(0.5d/ao)]2  km 


(222) 


where  h,  9  are  from  (21),  a  s  from  (15),  and  the  great-circle  path  distance,  d, 

\  i  c.  J 

is  a  starting  parameter  except  for  LOS  paths  where  it  is  calculated  from  (118), 


rLl  ,2  ^^,2  \l,2^+4^hl,2+ao^hLl,2+ao^sin^0,5dLl,2/ao^Z  km  ^223^ 
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where  hjj  2  are  from  (38,48) 


rBH  =  rLl  *  rL2  *  ds  kn' 


where  d$  is  from  (179) 


lr0  or  rWH  for  L0S 


r  =  greater  of(  or 


km 


fd  or  rBH  otherwise! 


where  rQ  is  from  (112),  and 


Lbf  =  32.45  +  20  log  (fr)  dB 


where  frequency,  f[MHz],  is  a  starting  parameter  [20,  p.  82], 


(2; 


(22? 


(226) 


9.  ATMOSPHERIC  ABSORPTION 

The  formulation  used  to  estimate  median  values  for  atmospheric  absorption  is 
similar  to  the  IF-73  method  [14,  Sec.  A. 4. 5].  Allowances  are  made  for  absorption 
due  to  oxygen  and  water  vapor  by  using  surface  absorption  rates  and  effective  ray 
lengths  where  these  ray  lengths  are  lengths  contained  within  atmospheric  layers 
with  appropriate  effective  thicknesses.  Geometry  associated  with  this  formulation 
is  shown  in  Figure  12  along  with  key  equations  relating  geometric  parameters.  This 
geometry  is  used  to  calculate  effective  ray  lengths  applicable  to  the  oxygen,  rgo, 
rain  storm,  r  for  Section  10.4,  and  water  vapor,  r  ,  layers  for  different  path 

Cj  CW 

configurations.  The  rain  storm  is  assumed  to  occur  between  the  facility  and  its 
maximum  LOS  distance,  d^  from  (43),  so  that  only  the  facility  horizon  ray  is  con¬ 
sidered  in  the  calculation  of  r  for  beyond- the- horizon  paths. 

i  es 

For  1 ine-of-sight  paths,  (d  £  d^ )  where  d  is  a  specified  parameter  except  in 
the  LOS  region  where  it  is  calculated  using  (118),  d^  is  from  (43),  the  Figure  12 
expressions  are  used  to  calculate  effective  ray  lengths,  r  ,  with  H  ,  «  =  H,  « 

GO ) S  jW  y  I  >  c  I  |  c 

from  (106),  for  earth,  a^,  =  a  from  (104),  and  p  =  eh1  from  (115a). 


Parameter  values  for  Hyl  km,  km,  and  a^  km  and  0  are  defined  In 
the  text  for  llne-of-slght,  single-horizon,  and  two-horizon  paths. 


Note:  Values  of  Te0jW  for  oxygen  and  water  vapor  are  taken  as  3.25  and 
1.36  km  [12,  Table  A. 2],  respectively.  The  value  of  T0S  used  to  estimate 
the  in-storm  ray  length,  Res,  for  rain  attenuation  is  dependent  on  the 
storm  site  (5ect1on  10.4). 


FIGURE  12.  GEOMETRY  ASSOCIATED  WITH  ATMOSPHERIC  ABSORPTION  CALCULATIONS. 


For  single  horizon  paths,  (d^  <  d  _<  d^  +  )  where  djj  is  from  Figure  5 

and  d^Q-j  is  from  (52),  the  Figure  12  expressions  are  used  with  one  (Rgs)  or  two 
(R„_  sets  of  starting  parameters  and  the  R  's  obtained  with  these  are 

GOjW  CU ,5jW 

caned  rleo,s,w  and  r2eo,W  In  the  calculation  of  r]eo  s  w,  Hy)  =  lesser  of 
(hei  or  hLrl)  and  Hy2  =  greater  of  (hgl  or  hj  rl )  where  hfil  is  from  (24)  and  hLrl  is 
from  (39);  a^  =  a  from  (16)  and  3  =  6^  from  (40)  if  =  hgl ,  otherwise  3  =  0|_ 
from  (42).  In  calculation  for  r2eQ  w>  H  j  =  lesser  of  (hg2  from  (24)  or  hLrl )  and 
H^2  =  greater  of  (hg2  or  hLr2)  with  3  =  9g2  from  (4?)  if  =  he2,  otherwise 
3  =  -9,  .  Values  for  r  are  then  obtained  using: 

l  QA  C  lit 


res  =  rles 


kin 


reo,w  rleo,w  +  r2eo,w 


(227a) 

(227b) 


For  two-horizon  paths  (dL]  +  dLQ1  <  d),  the  Figure  12  expressions  are  also  used 

with  one  (R  )  or  two  (R  )  sets  of  input  parameters.  The  results  obtained  are 
0  S  00  )  w 

called  rleo>s>w  and  r2e0>SjW>  where  (227)  is  used  to  determine  reQjS>w  values.  In 
calculations  for  r]e0  s  w,  H  j  =  lesser  of  {hel  or  hy  from  (201)),  =  greater 

of  (hgl  or  hv>,  3  =  0gl  if  H  ,  =  hgl ,  otherwise  3  =  -Tan_19A1  where  0A1  is  from  (196). 
In  calculations  for  r2eQ  s  w>  H ^  =  lesser  of  (hy  or  he2> ,  Hy2  =  greater  of  {hy  or 
he2),  a^  «  a  and  3  =  0e2  if  H  j  =  he2,  otherwise  3  =  -Tan-19A2  where  0A2  is  from  (196). 

Surface  absorption  rates  for  oxygen  and  water  vapor,  y00  w  are  used  with 

effective  ray  lengths,  reQ  to  obtain  an  estimate  for  atmospheric  absorption, 

Aa;  i.e.,: 


Y 


oo  eo 


+  Y, 


ow  ew 


dB 


(223) 


Values  for  y  are  obtained  by  interpolating  between  values  taken  from  the 
00  ,w 

Rice  et  al .  curves  [34,  p.  3-7]. 
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10.  VARIABILITY 


Total  variability,  Y^(q),  for  (1)  is  calculated  from: 

VE(q)  =  i  yv|(q)  ♦  rj(q)  +  Y*(q)  +  Y4(q)  dB  (229) 

+  for  q  £  0.5 
-  otherwise 


where  Ye(q)  is  the  variability  associated  with  long-term  (hour-to-hour)  power  fading 
(5ec.  10.1),  Y^q)  is  short-term  (within-the-hour)  fading  associated  with  surface 
reflection  multipath  (Sec.  10.2)  and  tropospheric  multipath  (Sec.  10.3),  Yr(q)  is 
rain  attenuation  fading  (Sec.  10.4),  and  Yj(q)  is  ionospheric  scintillations 
(Sec.  10.5).  When  the  long-term  (hourly  median)  time  availability  option  [20,  p.  103] 
is  selected,  short-term  fading  is  neglected  (i.e.,  Y7T(q )  =  Yr(q)  =  Yj (q )  =  0)  and 
the  resulting  variability  may  be  regarded  as  the  variability  of  hourly  medians. 

Since  Yit(0.5)  =  Yr(0.5)  *  Yj (0. 5)  *  0,  Y£(0. 5)  does  not  change  when  short-term 
fading  is  neglected.  For  q  <  0.98,  Yf(q)  =  0;  i.e.,  rain  attenuation  only  occurs 
during  2  percent  of  the  time.  Furthermore,  Y j (q )  *  0  except  for  earth/satellite 
paths  that  pass  through  the  ionosphere  which  is  at  an  altitude  of  about  350  km. 


10.1*  Long-Term  Power  Fading,  Yg(q) 

The  variability  associated  with  long-term  (hourly  median)  power  fading  for  (229) 
is  designated  Ye(q )  where  q  is  the  time  avai lability  parameter  of  Section  2.1  and 
the  sign  associated  with  Yfi(q)  values  is  such  th*t  the  positive  values  associated 
with  q  _<  0. 5  will  decrease  transmission  loss  or  increase  received  power  1 evel s .  All 
long-term  variability  options  [20,  p.  103]  use  an  effective  distance,  dg  [34,  Secs. 
10,  III -6 ,  1 1 1-7]  which  is  calculated  as  follows: 

dqs  =  65  (100/f )1/3  km  (230) 

where  f[MHz]  is  frequency: 


d  =  d.  +  d 
q  Lq  qs 


km 


(231) 
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where  d^  is  a  total  .smooth  earth  horizon  distance  determined  by  ray  tracing  (Sec.  3) 
with  N .  =  329  in  (17)  which  would  correspond  to  a  9000  km  effective  earth  radius 
[34,  p.  4-4] 


de 


130d/d  for  d  <  d 

q  -  q 


km 


(232) 


130  +  d  -  d  otherwise 
9 


where  d  is  the  gn  ic-circle  path  distance  and  is  a  specified  parameter  except  in  the 
LOS  region  e  it  is  calculated  via  (118).  Key  parameters,  g ( 0. 1  or  0.9,f),  V(0.5) 


and  Yq(0.1  or  0.9)  for  the  long-term  variability  normally  used  are  determined  as 


follows: 


0.21  sin[5. 22  log(f/200)]+l .28  tor  60<f<1600  MHz 


9(0.1, f)  = 


(233a) 


1.05  for  f  >  1600  MHz 


0.18  sin[5. 22  1 og ( f / 200 ) ]+l . 23  for  60<f<1600  MHz 


9(0.9, f)  = 


(233b) 


1.05  for  f  >  1600  MHz 


(fm  -  f„)  exp(-c2d^2)  (234) 


clde'  "  f  2  ]  exp(‘c3d"3)  +  f2  dB  (235) 

where  the  values  used  for  the  parameters  C|,  c^,  c^,  n^ ,  n^,,  n3,  fm,  and  f^  depend 
on  whether  V(0.5)  [34,  Table  III. 5,  Climate  1],  Y(0.1)  [34,  Table  III. 3,  all  hours 
all  year],  or  Y(0.9)  [34,  Table  III. 4,  all  hours  all  year]  is  calculated.  This 
selection  is  based  on  a  recommended  model  [11,  p.  19]  that  was  tested  against 
air/ground  data  [10,  Sec.  4.3].  However,  other  options  such  as  different  time 
blocks  (Sec.  10.1.1),  climates  (Sec.  10.1.2),  or  a  mix  performed  to  meet  particular 


fn  = 


and 


V(0. 5) 

¥0(°.l) 


A 


-Yo 


(0.9) 


} 


conditions  (Sec.  10.6)  may  be  used.  The  key  parameters  from  (233,235)  or  elsewhere 
(Secs.  10.1.1,  10.1.2)  are  used  to  obtain  the  final  long-term  variability  distribu¬ 
tion  as  follows: 


Y(0.1  or  0.9)  *  g(q  =  0.1  or  0.9,f)  YQ(q  =  0.1  or  0.9)  dB 


C cY(0. 1 )  for  q  <  0. 5  ^ 
Y(q)  *  <  0  for  q  = 

^c.Y(0.9)  for  q 


1  -  0.5  > 

>  0.5  J 


dB 


where  c  values  for  the  q  desired  are  selected  from  Tables  [15,  p.  34], 

fo.5  -  tt"1  Tan-1  [20  log  (328hl)]  for  LOS  *1 

f9h  *  *  paths  with  >  0  r 

J  otherwise  *  J 


(236) 


(227) 


(238) 


where  ohl  is  from  (115)[13,  p.  8;  34,  p.  III. 43]. 


Vei(q)  -  f0hKq> 


dB 


Ve(0*5)  =  W(0*5)  dB 


Yt  =  Lb(0.5)  -  [Lbf  -  Vfi(0.5)  -  20  log  (RTg  +  RTc)] 
where  Lb(0. 5)' is  from  (2),  is  from  (226),  R-j.g  c  are  from  (119), 


dB 


AYI  =  ^ 


0  if  lobing  option  [20,  p.  99]  is  used 
and  aircraft  is  within  10  lobes  of 
its  radio  horizon 

(Lbf  -  3)  -  CLbr  -  YeI(0-U]  otherwise 


> 


(239) 

(240) 

(241) 


(242) 


where  Lbr  is  from  (3), 
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LV- 


$ 

’  v 

rV 

1 


9 


r:-: 

I 


I 


$ 


0  if  Ayj  £  0 

Ay  =  ^  10  if  Ayj  >_  10 


(243) 


YI 


otherwise 


r 


Ye(q<0.1)  =  < 


lesser  of  CYg j ( q )  or  Y^]  for  lobing 
lesser  of  [Yg j (q )  or 


r  dB 


L,  +  fi.  - 
br 

otherwise 


uf  ”  CY^ 


(244) 


J 


where  c^  is  6,  5.8,  and  5  dB  for  q  values  of  0.0001,  0.001,  and  0.01,  respectively. 


Ye(q>0.1)  =  YeI(q)  dB  . 


(245) 


10.1.1  Time  Blocks 

Long-term  variability  options  for  IF-77  include  variabilities  appropriate  for 
the  time  blocks  shown  in  Table  1  [20,  p.  103].  These  blocks  and  seasonal  groupings 
are  used  to  describe  the  diurnal  and  seasonal  variations  in  a  continental  temperate 
climate  [34,  Sec.  III. 7.1].  They  are  incorporated  into  (235)  by  selecting  appro¬ 
priate  constants  from  Rice  et  al .  [34,  Tables  III. 2,  III. 3,  and  III. 4].  The 
expressions  for  g(0. l,f)  and  g(0.9,f)  given  in  (233)  are  used  for  time  block  varia¬ 
bilities. 

If  a  combination  of  time  blocks  is  appropriate,  various  distributions  can  be 
mixed  (Sec.  10.6). 

10.1.2  Climates 

Options  to  use  various  climates  are  included  in  IF-77  [20,  p.  103];  i.e., 

(1)  equatorial,  (2)  continental  sub-tropical,  (3)  maritime  sub-tropical,  (4)  desert, 
(5)  mediterranean,  (6)  continental  temperate,  (7a)  maritime  temperate  overland, 

(7b)  maritime  temperate  oversea,  and  (8)  polar.  The  formulation  used  is  based  on 
algebraic  expressions  fitted  to  modified  versions  of  curves  provided  in  CCIR 
Report  238-4  [7]  by  Hufford  and  Longley  [D0C-BL,  informal  communication; 

15,  Sec.  4.3;  29,  Sec.  4.4.25]  and  may  be  summarized  as  follows: 
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TABLE  1.  TIME  BLOCK  RANGES  [34,  Sec.  1 1 1. 7.1] 


NO. _  MONTHS  HOURS 


1 

Nov. 

-  Apr. 

0600  -  1300 

2 

Nov. 

-  Apr. 

1300  -  1800 

3 

Nov. 

-  Apr. 

1800  -  2400 

4 

May 

-  Oct. 

0600  -  1300 

5 

May 

-  Oct. 

1300  -  1800 

6 

May 

-  Oct. 

1800  -  2400 

7 

May 

-  Oct. 

0000  -  0600 

8 

Nov. 

-  Apr. 

0000  -  0600 

Summer 

May 

-  Oct. 

ALL-HOURS 

Winter 

Nov. 

-  Apr. 

ALL-HOURS 

V ( 0 . 5) 
Yq(0.1) 
Yo(0.9)- 


<de/bl>‘ 


1  +  (dg/bi V 


C1  + 


1  +  [(de  -  b2)/b3]t 


dB 


where  dg  is  from  (232),  and  appropriate  values  for  the  b1 ,  b^,  b3,  c1 
used  [15,  p.  33], 


g(o.l.f) 


r 


1  for  all  climates  except  2,  4,  and  6, 
0.18  sin  [5  log  (f/200)]  +  1.06 


< 


for  60  _<  f  £  1500  MHz  in  Climates  2  and  6, 
1  suggested  for  60  <_  f  200  MHz  in  Climate  4, 
0.10  sin  [5  log  (f/200)]  +  1.02, 


for  200  ^  f  <  1500  MHz  in  Climate  4, 

0.93  for  f  >  1500  MHz  In  Climates  2,  4,  and  6 


1  for  all  climates  except  6, 
0.13  sin  [5  log  (f/200)]  +  1.04 


9(0.9,f) 


< 


for  60  f  £  1500  MHz  in  Climate  6,  V 


0.92  for  f  >  1500  MHz  in  Climate  6 


(247b) 


where  f[MHz]  is  an  IF-77  input  parameter  [20,  p.-  82].  When  one  of  these  climate 
options  is  selected,  the  key  parameters  obtained  from  (246,247)  are  used  in  (236,240). 
If  a  mixture  of  two  or  more  climates  is  appropriate,  various  distributions  can  be 
mixed  (Sec.  10.6). 


10.2  Surface  Reflection  Multipath 

Multipath  associated  with  reflections  from  the  earth's  surface  is  considered 
as  part  of  the  short-term  (within  the  hour)  variability  for  1 ine-of-sight  paths, 
and  is  used  only  when. the  time  availability  option  for  "instantaneous  levels 
exceeded"  is  selected  [20,  p.  103].  Contributions  associated  with  both  specular  and 
diffuse  reflection  components  may  be  included  though  the  specular  component  is  not 
allowed  to  make  a  full  contribution  when  it  is  also  used  in  determining  the  median 
levels  (e.g.,  when  lobing  option  is  selected  [20,  p.  S?]).  These  contributions  are 
incorporated  into  the  variability  part  of  the  model  via  the  relative  power  level, 

WR,  which  is  used  to  compute  Y^q)  for  (229)  as  shown  in  Section  10.3.  Formulas 
used  to  calculate  WB  may  be  summarized  as  follows: 


1  if  Ay  <  0 


AY 


0.1  if  Ay  >  9 


0.5  [1.1  +  0.9  cos  (iTAy/9 ) ]  otherwise 


(243) 


where  Ay  is  from  (243), 


10.3  Tropospheric  Multipath 

Tropospheric  multipath  is  caused  by  reflections  from  atmospheric  sheets  or 
elevated  layers,  or  additional  direct  (nonref lected)  wave  paths  [2;  9,  Sec.  3.1], 
and  may  be  present  whe:i  antenna  directivity  is  sufficient  to  make  surface  reflections 
negligible.  It  is  considered  as  part  of  the  short-term  (within  the  hour)  variabil- 
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ity  used  only  when  the  time  availability  option  for  "instantaneous  levels  exceeded" 
is  selected  [20,  p.  103],  and  incorporated  into  the  variability  part  of  the  model 
via  the  relative  power  level,  Wfl.  This  is  used  to  compute  Y  ( q )  for  (229)  as  shown 
in  the  latter  part  of  this  section.  The  formulation  for  Wa  may  be  summarized  as 

a 

follows: 


10  log  (f 


-i>  - 


84.26  for 


d  —  dML 


F  = 


and  is  not  calculated  otherwise 


I 


dB 


(253) 


where  frequency,  f[MHz]  is  an  input  parameter  [20,  p.  82],  rgw  is  calculated  as 
indicated  on  Figure  12,  d  is  from  (118)  for  d  <_  d^ ,  and  d^  is  from  (43). 


10S 


40  dB  for  F  <  0.14 
-20  dB  for  F  >18.4 
otherwise  obtained  from  curves 


dB 


(254) 


where  curves  for  the  Nakagami-Rice  probability  distribution  [34,  p.  V-8]  are  used  by 
selecting  the  "K"  (which  becomes  K^g^)  that  corresponds  to  a  "Yif(0.99)"  of  -F 


MKa  =  (-20-KML)/0. 02618  dB/rad 


(255) 


where  is  evaluated  at  the  maximum  LOS  range  (i.e.,  d  =  d^). 


K 


t 


kL0S  for  d  1  dML 

-20  for  9S  >  0.02618  rad 

Kj^L  +  ^Kaes  otherwise 


dB 


(256) 


where  is  the  scattering  angle  from  (180)  and  is  negative  in  the  LOS  region. 


-K./10 

W  *  10  1 

a 


(257) 
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Relative  powers  and  WR  from  (252)  are  combined  to  determine  K,  which  is  the 
ratio  in  decibels  between  the  steady  component  (e.g.,  direct  ray),  arid  the  Rayleigh 
fading  component  (e.g.,  surface  reflection  and  tropospheric  multipath)  using: 

K  -•  -10  log  (WR  +  Wa)  dB  (256) 

The  Nakagami-Rice  probability  distribution  for  Y^q)  of  (229)  is  then  selected  from 
the  Rice  et  al.  curves  [34,  p.  V-8]  by  using  K. 

10.4  Rain  Attenuation 

The  rain  attenuation  model  used  in  IF-77  is  largely  based  on  material  in 
informal  papers  by  C.  A.  Samson  (DOC-BL)  on  "Radio  Propagation  Through  Precipitation" 
and  "Rain  Rate  Distribution  Curves."  This  discussion  is  a  shortened  version  of  the 
description  previously  provided  [15,  Sec.  4.4],  and  the  maps  and  tables  provided 
there  are  not  repeated  here. 

Two  options  for  rain  attenuation  are  available  in  IF-77.  The  first  is  for  use 
in  a  "worst  case"  type  analysis  where  a  particular  rainfall  attenuation  rate  is 
assumed  for  the  in-storm  path  length,  and  the  additional  path  attenuation  associ¬ 
ated  with  rain  is  simply  taken  as  the  product  of  this  attenuation  rate  (in  dB/km) 
and  the  in-storm  ray  length  [20,  p .  94] .  This  ray  length  is  determined  in  accord¬ 
ance  with  the  method  discussed  in  step  4  of  option  two. 

Option  two  involves  computer  input  of  rain  zone  (which  determines  a  rainfall 
rate  distribution)  and  storm  size  [20,  p.  94].  Storm  size  (diameter  or  long 
dimension)  is  assumed  to  be  one  of  three  options:  5,  10,  or  20  km  (corresponding 
approximately  to  a  relatively  smail,  average,  or  very  large  thunderstorm).  The 
maximum  distance  used  in  calculating  path  attenuation  with  this  option  is  the  storm 
size  since  it  is  assumed  that  only  one  storm  is  on  the  path  at  a  time.  The  process 
used  to  include  rain  attenuation  estimates  i n; IF-77  for  this  option  may  be  summarized 
as  follows: 

1.  Determine  point  rain  rates.  Point  rain  rates  (rate  at  a  particular  point 
of  observation)  not  exceeded  for  specific  fractions  of  the  time  are  determined  for 
the  rain  zone  of  interest  [15,  p.  38]. 

2.  I)etermine_path  average  rain  rates.  Each  point  rain  rate  resulting  from 
step  1  is  converted  to  a  path  average  rain  rate  by  using  linear  interpolation  to 
obtain  a  multiplying  factor  [15,  p.  39]. 
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3.  Determine  attenuation  rate.  For  each  path  average  rain  rate  resulting 
from  step  2,  an  attenuation  rate  Arr(q)[dB/km]  is  determined  using  linear  interpola¬ 
tion  between  the  values  previously  determined  for  various  rainfall  rates  and  frequency 
[15,  p.  40]. 

4.  Determine  the  in-storm  ray  length.  First,  the  length  of  the  direct  ray 
rgs  that  is  within  Tes  of  the  earth's  surface  is  determined  from  (227)  by  using  the 
method  described  in  Section  9  where  Te$  is  taken  as  the  storm  size;  i.e.,  storm 
height  is  assumed  to  be  equivalent  to  storm  diameter.  Then,  the  final  in-storm  ray 
length,  r$,  is  calculated  using: 


rs  = 


f  ^es  ^  res  -  ^es 


■\ 


r  otherwise 
es 


(259) 


a 


5.  Determine  rain  attenuation  values.  Values 
particular  fraction  of  time  are  calculated  using 


for  the  attenuation,  A^(q)  for 


A 

r 


(q)  = 


0  for  q  _<  0. 98 
Arr(q)rs  otherwise 


dB 


(260) 


where  Ar,  ^q)  values  come  from  step  3  and  the  value  for  r$  is  from  step  4.  Distribu¬ 
tions  of  rain  attenuation  are  zero  for  q  <  0.98  [js,  p.  38]. 

6.  Combine  rain  attenuation  variability  with  other  variabilities.  Variability 
for  rain  attenuation,  Yr(q),  is  related  to  the  distribution  of  rain  attenuation  by: 

Yr(q)  =  -Ar(q)  (261) 


It  is  combined  with  other  variabilities  in  (229). 
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10.5  Ionospheric  Scintillation 

Variability  associated  with  ionospheric  scintillation,  Yj(q)  for  (229),  for 
paths  that  pass  through  the  ionosphere  (i.e. ,  on  earth/satellite  paths)  at  an  alti¬ 
tude  of  about  350  km  is  included  in  IF-77.  This  variability  may  be  specified 
directly  by  the  selection  of  a  scintillation  index  group  or  by  using  a  weighted 
mixture  of  distributions  where  the  weighting  factors  are  estimated  for  specific 
problems  [20,  p.  91;  15,  Sec.  4.5].  Provisions  are  included  to  allow  Yj(q)  to 
change  with  earth  facility  latitude  when  a  geostationary  satellite  is  involved  and 
the  earth  facility  locations  are  alcng  the  subsatellite  meridian.  When  this  pro¬ 
vision  is  used,  Y-j (q )  is  obtained  from  previously  prepared  136  MHz  data 
[15,  p.  45]  and  Yj(q)  is  calculated  as  follows: 

r 


n  = 


a 


for  8pL  <  17° 


or  epL  >  52° 


(eFL  -  17)/7 


for 


17°  <  en  < 

rL 


24° 


> 


<  0pL  <  45° 


for  24° 

-  (52  -  0pL )/ 7  for  45°  <  epL  <  52c 


(262) 


J 


where  6pL  is  the  magnitude  of  the  earth  facilities  latitude  in  degrees,  and 


Yj(q)  =  (136/f)nY136(q)  dB 


(263) 


r  though  this  scaling  factor  is  built  into  the  programs  that  utilize  IF-77,  only 
minor  program  modifications  would  be  required  to  use  other  simple  scaling  methods. 

In  addition,  the  distribution  mixing  methods  of  'Section  10.6  could  be  used  to  create 
Yj(q)'s  applicable  to  specific  situations. 

10.6  Mixing  Distributions 

Subroutines  have  been  incorporated  into  the  IF-77  computer  programs  to  allow 
the  distributions  that  characterize  portions  of  the  variability  associated  with  a 
particular  model  component  to  be  mixed  in  order  to  obtain  the  total  variability  for 
that  component.  For  example,  different  fractions  of  the  time  may  be  characterized 


by  signal  level  distributions  associated  with  different  ionospheric  scintillation 
groups;  and,  with  these  subroutines,  they  can  be  weighted  and  combined  (mixed)  to 
obtain  the  total  variability  associated  with  ionospheric  scintillations  (Sec.  10.5). 

The  process  of  mixing  N  cumulative  variability  distributions  may  be  summarized 
as  follows; 

1.  Select  M  (ten  or  more)  levels  of  variability  V-| ,  ...,  ,  ...,  VM  that 

cover  the  entire  range  of  the  transmission  loss  (or  power  available,  etc.)  values 
involved. 

2.  Determine  the  fraction  of  time  (weighting  factor)  for  which  each  distri¬ 
bution  is  applicable;  i.e.,  W-j ,  ...,  Wj,  ...,  W^. 

3.  Determine  the  tire  availability  (fraction  of  time  during  which  a  distri¬ 

bution  is  applicable  that  a  specific  level  of  transmission  loss  is  not  exceeded) 
for  each  distribution  at  the  selected  levels;  i.e.  ,<  •  •  • »  . . . »  q^ . 

4.  Calculate  time  availabilities  for  the  mixed  distribution  that  corresponds 
to  the  variability  levels  selected,  i.e.; 

ql  =  qllWl  +  **•  +  qljWj  +  •**  +  qlNWN 


qi  =  qilWl 


+  +  qijWj  + 


+  qiNWN 


(264) 


qM  "  qMlW1  +  +  qMjWj  +  +  qMNWN 

This  process  is  the  same  as  the  one  used  by  Rice  et  al .  [34,  Sec.  111.7.2]  to 
combine  transmission  loss  distributions  for  time  blocks  to  obtain  distributions  for 
summer  and  winter.  It  is,  also,  essentially  the  same  as  the  method  recommended  by 
Whitney  et  al.  [38,  p.  1099;  39,  Sec.  6]  to  combine  distributions  of  fading  associ¬ 
ated  with  various  ionospheric  scintillation  index  g  oups  (Sec.  10.5). 


■*! 
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When  this  process  is  used  to  mix  distributions  of  long-term  variability,  the 
required  variability  functions  are  obtained  from: 

Vc(q)  =  [V(0.5)  +  Y(q)]|c  dB  (265) 

where  . ..|c  indicates  that  the  V( 0.5)  and  Y(q)  are  appropriate  for  the  conditions 
(e.g.,  time  block  or  climate)  associated  with  a  particular  value  of  the  subscript  c. 
For  example,  V ( 0. 5 )  and  Y (q )  values  for  different  climates  can  be  obtained  by  using 
247)  with  (236,  237)  and  mixing  can  be  used  to  estimate  variability  for  areas  near  a 
border  between  two  different  climate  types.  After  mixing,  Y(q)  values  needed  for 
(239)  may  be  obtained  by  using: 

Y(q)  =  V ( q )  -  V(0. 5)  dB  (266) 

where  all  variables  in  (266)  are  associated  with  the  resulting  mixed  distribution. 
Similarly,  when  mixing  variabilities  associated  with  ionospheric  scintillation, 

VIc(q)  -  V,(q)|c  dB  (267) 

and  the  distribution  resulting  from  the  mixing  is  taken  as  Y j ( q )  for  later  calcu¬ 
lations. 


1 1 .  SUMMARY 

The  IF-77  electromagnetic  wave  propagation  model  was  discussed,  and  references 
were  provided  so  that  more  information  on  r.pecific  items  could  be  obtained.  A  brief 
description  of  the  model  provided  in  Section  2  is  followed  by  a  systematic  discus¬ 
sion  of  model  components.  Readers  with  a  further  interest  in  IF-77  are  encouraged 
to  obtain  a  copy  of  the  "Applications  Guide  for  Propagation  and  Interference 
Analysis  Computer  Programs  (0.1  to  20  GHz)"  [20]. 
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APPENDIX 


ABBREVIATIONS,  ACRONYMS,  AND  SYMBOLS 

This  list  includes  most  of  the  abbreviations,  acronyms,  and  symbols  used  in 
this  report.  Many  are  similar  to  those  previously  used  [14,  15,  20,  34].  The 
units  given  for  symbols  in  this  list  are  those  required  by  or  resulting  from  equa¬ 
tions  as  given  in  this  report  and  are  applicable  except  when  other  units  are 
specified. 

In  the  following  list,  the  English  alphabet  precedes  the  Greek  alphabet,  letters 
precede  numbers,  and  lower-case  letters  precede  upper-case  letters.  Miscellaneous 
symbols  and  notations  are  given  after  the  alphabetical  items. 

a  Effective  earth  radius  [km]  calculated  from  (16). 

a  An  adjusted  effective  earth  radius  [km]  shown  in  Figure  8, 

from  (104). 

aQ  Actual  earth  radius,  6370  km,  to  about  three  signficant  figures. 

a  An  effective  earth  radius  [km]  used  in  Figure  12,  and  defined 

Y  for  different  path  types  in  Section  9. 

a1  ^  Effective  earth  radii  [km],  from  (62). 

a^  4  Effective  earth  radii  [km],  from  (63). 

A  A  parameter  used  in  tropospheric  scatter  calculations, 

from  (210). 

Ag  Atmospheric  absorption  [dB],  from  (228). 

A A  Combined  diffraction  attenuation  [dB]  at  d  =  d^,  from  (97). 

A^  Attenuation  [dB]  associated  with  diffraction  beyond  the 

horizon,  from  (100). 

A^q  Intercept  [dB]  for  the  beyond-the-horizon  combined  diffraction 

1  attenuation  line,  from  (99). 

Adx  Ad  [dB]  at  dx,  discussed  after  (221). 

A  .  Angles  [rad]  defined  and  used  in  Figure  12  only. 

6  |  L 
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Attenuation  line  intercept  [dB]  of  rounded  earth  diffraction  for 
path  F-O-ML  (Figure  7),  from  (54,78). 

Effective  area  [dB-sq  in]  of  an  isotropic  antenna,  from  (9). 

Knife-edge  diffraction  attenuation  [dB]  for  path  F-O-A 
(Figure  7) ,  from  (93) . 

Knife-edge  diffraction  attenuation  [dB]  for  path  F-O-ML 
(Figure  7),  from  (87). 

Knife-edge  diffraction  attenuation  [dB]  at  d^,  from  (168). 

Line-of-sight  attenuation  [dB],  from  (177). 

A  parameter  [km]  used  in  tropospheric  scatter  calculations, 
from  (183) . 

Combined  diffraction  attenuation  [dB]  at  d^L,  from  (95). 

ARq  [dB]  at  dQ,  from  (175). 

Attenuation  line  intercept  [dB]  of  rounded  earth  diffraction 
path  0-A  (Figure  7),  from  (54s78). 

Rounded  earth  diffraction  attenuation  [dB]  for  path  p, 
from  (54,78). 

Arp  [dB]  at  d^,  from  (85). 

Rounded  earth  diffraction  attenuation  [dB]  for  path  F-O-A 
(Figure  7),  from  (85). 

Arp  [dB]  at  d^L,  from  (85). 

Rounded  earth  diffraction  attenuation  [dB]  for  path  0-A 
(Figure  7)  at  distance  d^,  from  (86). 

Attenuation  [dB]  in  the  line-of-sight  region  where  the  diffrac¬ 
tion  effects  associated  with  terrain  are  negligible  (d  <  d  ), 
from  (175).  0 

Attenuation  [dB]  due  to  rain  for  a  fraction  of  time  q, 
from  (260). 

Attenuation  rate  [dB/km]  associated  with  rain  and  a  fraction  of 
time  q,  (Sec.  10.4,  Step  3). 

Arp  [dB]  at  d5,  from  (169). 

Forward  scatter  attenuation  [dB],  from  (220). 


A  [dB]  at  d  ,  discussed  after  (221). 

5  A 

Terrain  attenuation  [dB],  from  (221). 

A  conditional  aajustment  factor  [dB]  used  to  prevent  available 
signal  powers  from  exceeding  levels  expected  for  free-space 
propagation  by  unrealistic  amounts,  from  (243). 

An  initial  value  of  Ay  [dB],  from  (242). 

Rounded  earth  diffraction  attenuations  [dB],  from  (76). 

Combined  diffraction  attenuation  [dB],  at  d^,  from  (170). 

Long-term  variability  parameters  with  values  previously 
provided  [15,  p.  33],  used  in  (246). 

A  parameter  used  in  the  G_  formulation  of  (83). 

hi, 2 

Parameters  used  in  rounded  earth  diffraction  formulation, 
from  (79). 

A  parameter  used  in  forward  scatter  formulation,  from  (217). 


Parameters  used  in  the  rounded  earth  diffraction  formulation, 
from  (68). 

A  parameter  used  in  (237)  with  values  taken  from  tables 
[15,  p.  34]. 

Phase  of  plane  earth  reflection  coefficient  relative  to  tt  for 
eliptical,  horizontal,  and  vertical  polarization. 

A  variability  limiting  parameter  [dB]  used  in  (244)  and  defined 
just  after  (244). 

Long-term  variability  parameters  for  (234,235,246)  that  are 
discussed  following  (235). 

An  exponential  atmosphere  parameter,  from  (18). 


A  parameter  used  in  tropospheric  scatter  calculations, 
from  (218). 

A  Fresnel  integral  [34,  1 1 1-18] ,  for  (92,137,167). 

Great-circle  distance  between  facility  and  aircraft.  For 
1 ine-of-sight  paths,  it  is  calculated  via  (113). 

Decibel,  10  log  (dimensionless  ratio  of  powers). 

Antenna  gain  in  decibels  greater  than  isotropic. 


dB/km  Attenuation  [dB]  per  unit  length  [km]. 

dB-sq  m  Units  for  effective  area  in  terms  of  decibels  greater  than 

an  effective  area  of  1  square  meter;  i.e.,  10  log  (area  in 
square  meters). 

dBW  Power  in  decibels  greater  than  1  watt. 

dB-W/sq  m  Units  of  power  density  in  terms  of  decibels  greater  than  1  watt 

per  square  meter;  i.e.,  10  log  (power  density  expressed  in  watts 
per  square  meter). 

dN  A  parameter  [per  km]  used  in  tropospheric  scatter  calculations, 

,  from  (184). 

d^  A  facility-to-aircraft  (Figure  7)  distance  [km],  from  (96). 

d  Counterpoise  diameter  [km],  a  model  input  parameter  [20,  p.  88, 

FACILITY  ANTENNA  COUNTERPOISE  DIAMETER]. 

dd  Initial  estimate  of  dQ,  from  (171). 

d  Effective  distance  [km],  from  (232). 


JLE1 


LM 


JL01 ,2 


JLp 


Lpl  ,2 


Lq 


A  distance  [km]  used  in  facility  horizon  determination, 
from  (35). 

An  initial  value  for  the  facility  horizon  distance  [km]  that  is 
based  on  affective  earth  radius  geometry,  and  shown  in  Figure  4. 
It  may  be  specified  [20,  p.  90]  or  calculated  as  indicated  in 
Figure  5,  f>om  (33,36). 

Maximum  distance  [km]  for  which  the  facility-to-aircraft  path 
has  a  common  horizon,  from  (45). 

Smooth  earth  Horizon  distances  [km]  for  path  0-A  (Figure  7), 
from  (52,53). 

Total  horizon  distances  [km]  for  path  p,  from  (59). 

Radio  horizon  distances  [km]  for  path  p,  from  (52,53). 

Smooth  earth  horizon  distances  [km]  determined  via  ray  tracing 
(Sec.  3)  over  a  9000-km  (4860  n  mi)  earth,  discussed  after 
(231). 


LR2 


Horizon  distance  for  aircraft  shown  in  Figure  4  and  discussed 
preceding  (43). 


—4  mii  « 
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d^s  Total  smooth  earth  horizon  distance  [km],  from  (28). 

d.  .  The  sum  of  the  smooth  earth  distances  of  path  0-A  (Figure  7), 

LSM  from  (88). 

d.  n  o  A  smooth  earth  horizon  distance  [km]  for  the  facility  or 

aircraft  that  is  based  on  effective  earth  radius  geometry, 
from  (26). 

d.  R1  ?  Smooth  earth  horizon  distances  determined  via  ray  tracing 

’  (Sec.  3.1),  and  shown  in  Figure  3. 

dLsl  2  Smooth  earth  horizon  distances  determined,  from  (27). 

d..  Facility-to-horizon  distance  [km],  determined  as  shown  in 

Figure  5,  from  (41). 

d^2  Horizon  distance  for  aircraft,  from  (46). 

dL5  A  distance  [km],  calculated  from  (161). 

dMI  Maximum  1 ine-of-sight  distance  [km]  shown  in  Figure  4, 

ML  from  (43). 

d  The  largest  distance  [km]  in  the  1 ine-of~sight  region  at  which 

diffraction  effects  associated  with  terrain  are  considered 
negligible,  from  (172). 

dq  A  distance  [km],  calculated  from  (231). 

d  Distance  [km]  beyond  the  radio  horizon  at  which  diffraction  and 

q  scatter  attenuation  are  approximately  equal  for  a  smooth  earth, 

from  (230). 

ds  Distance  [km]  between  horizons,  from  (179). 

d  |  Smooth  earth  horizon  distance  [km]  for  the  obstacle  height  as 

shown  in  Figure  6,  from  (44). 

d  A  distance  [km]  just  beyond  the  radio  horizon  where 

As  £  20  dB  and  M&  £  M^,  discussed  after  (221). 

d  .  ?  Distances  [km]  used  in  tropospheric  scatter  calculations, 

zuc  from  (181,182). 


d71  „  Distances  [km]  used  in  tropospheric  scatter  calculations, 

’  from  (199,200). 


do  a  Distances  [km]  used  in  rounded  earth  diffraction  calculation, 

’  from  (60,61), 

d^  A  distance  [km]  calculated  from  (162). 

d^/g  The  largest  distance  [km]  at  which  a  free-space  value  of  basic 

transmission  loss  is  obtained  in  two-ray  model  of  reflection 
from  a  smooth  earth  with  an  effective  reflection  coefficient  of 
-1.  This  occurs  when  the  path  length  difference,  Ar  from  (114), 
is  equal  to  A/6. 

DOC-BL  United  States  department  of  Commerce,  Boulder  Laboratories. 

DOT  United  States  Department  of  Transportation. 

D/U(q)  Desired-to-undesired  signal  ratio  [dB]  exceeded  for  at  least  a 

fraction  q  of  the  time.  These  values  may  represent  instantaneous 
levels  or  hourly  median  levels  depending  upon  the  time  availa¬ 
bility  option  selected  [20,  p.  103]  and  are  calculated  via  (11). 

Dy  Divergence  factor,  from  (122). 

Dj  ^  Distances  [km]  shown  in  Figure  8  and  calculated  via  (109). 

e  2.718281828. 

eq.  Equation. 

exp(  )  Exponential;  e.g.,  exp(2)  =  e2  or  R  exp(-j<j>)  =  R  is  a 

phasor  with  magnitude  R  and  a  lag  of  <J>  radians.  e 

EIRP  Equivalent  isotropically  radiated  £ower  [dBW]  as  defined  by  (7). 

EIRPG  EIRP  [dBW]  increased  by  the  main  beam  gain  [dBi]  of  the 

receiving  antenna  as  in  (6). 

f  Frequency  [MHz],  an  input  parameter  [20,  p.  82], 

fss  Facil ity  £ite  surface  (Figure  2). 

f  ,  A  parameter  used  in  G  weighting  factor  and  calculated  from  (81). 

ci.z  h 

f  f-  Knife-edge  diffraction  loss  factors  determined  using  Fresnel 

9,C,V’D  integrals,  from  (137,92,167). 

f  v  Parameters  used  in  the  normally  used  variability  formulation 

and  discussed  following  (235). 

f^  Elevation  angle  correction  factor,  from  (238). 
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F 


FAA 


Fade  margin  [dB],  from  (?53). 
Federal  Aviation  Administration. 


A  specular  reflection  reduction  factor  associated  with  A 
from  (248). 


F.  .  Reflection  reduction  factor  ac'ociated  with  diffuse  reflection 

0  and  surface  roughness  [5, Figure  4],  from  (128). 

F^  A  factor  calculated  from  (173). 


F  Reflection  reduction  factor  associated  with  ray  lengths, 

from  (123). 

F^  2  Parameters  [dB]  calculated  from  (75). 


Fj  2  Parameters  [dB]  calculated  from  (74). 


F^r  Reflection  reduction  factor  associated  with  Ar,  from  (249). 

F  .  Specular  reflection  reduction  factor  associated  with  surface 

roughness,  from  (127). 


g(0. l,f)  or  Frequency  gain  factor,  from  (233,247). 
g(0.9,f) 


Voltage  gain  [V/V]  factors  associated  with  direct  and  reflected 
rays,  from  (141 ,142). 


gD1  2  Voltage  gain  [V/V]  of  terminal  antennas  in  the  direction  of  the 

’  direct  ray  (Figure  11),  relative  to  main  beam  gain. 


9hDl  ,2 


Voltage  gain  [V/V]  similar  to  g^,  2*  but  specifically  for 
horizontal  polarization.  * 


Voltage  gain  [V/V]  similar  to  gR1  but  specifically  for 
horizontal  polarization.  ’ 


Rg>c 


Gain  factors  used  in  (119)  for  earth  or  counterpoise  reflected 
rays,  from  (149). 


gRh 


Cain  factor  for  the  reflected  ray  and  horizontal  polarization, 
from  (144). 


Rv 


Gain  factor  for  the  reflected  ray  and  vertical  polarization, 
from  (143). 


g.,,  Voltage  gain  [V/V]  of  terminal  antennas  in  the  direction  of  the 

,',,i  reflected  ray  (Figure  11)  relative  to  main  beam  gain,  from 

(140). 


;Vvvv-;v  •’ 

f 

m 

9vdi  ,2 

Voltage  gain  [V/V]  similar  to  gD1  2»  but  specifically  for 
vertical  polarization.  'c 

M  9vR1,2 

Voltage  gain  [V/Vj  similar  to  gR,  2,  but  specifically  for 
vertical  polarization.  K  'c 

GHz 

Gigahertz  (109  Hz). 

get,r 

Gain  [dBi]  of  the  transmitting  or  receiving  antenna  at  an 
appropriate  elevation  angle,  from  (4). 

p  Ghl,2 

Values  [dB]  for  the  residual  height  gain  function,  from  (83).  I 

eI  6ff1,2 

G|ipl  2  ^or  Pat^  F-0-ML  (Figure  7),  from  (56,84). 

■v  Gh01,2 

r-'j 

G  [dB]  for  path  0-A  (Figure  7),  from  (56,84). 

Fpl  ,2 

•y  G 

y  hpl ,2 

Values  [dB]  for  the  weighted  residual  height  gain  function 
(Sec.  4.2)  for  path  p,  from  (56,84). 

W  gnt,r 

<  ^  r 
s.""n 

Normalized  gain  of  the  transmitting  or  receiving  antenna  at 
appropriate  elevation  angle  in  decibels  greater  tha j  the  main 
beam  (Gy  or  GR),  from  (145). 

i  Gr!-2 

<r 

Gains  [dB]  gR1  2  expressed  in  decibels,  from  (139). 

w 

GT,R 

Gain  [dBi]  of  the  transmitting  or  receiving  main  beam  for 
(4,6,7). 

n\ 

>  GW1 ,2 

Weighting  factor,  from  (82). 

R  G1 ,2,3,4 

Parameters  [dB]  calculated  from  (71). 

p  h 

Ray  height  [km]  above  msl  used  in  (17). 

i 

Height  [km]  of  facility  counterpoise  above  ground  at  the 
facility  site,  an  input  parameter  [20,  p.  88]. 

£  he 

An  effective  height  [km],  from  (34). 

c-'  W 

Effective  aircraft  altitude  [km]  above  msl,  from  (160). 

?  h 

hepl,2 

Effective  antenna  heights  [km]  for  path  p,  from  (50,51). 

:K  hel  ,2 

s 

Effective  antenna  height  [km]  above  hr  (Figure  3),  from  (24), 
and  shown  in  Figure  3, 

> 

•■v 

*■«' 

s* 
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hf  Height  [km]  of  the  facility  antenna  above  its  counterpoise 

L  (Figure  2). 

hLE1  An  initial  value  for  the  facility  horizon  elevation  [km]  above 

msl  that  is  based  on  effective  earth  radius  geometry  (Figure  4). 
It  may  be  specified  [20,  p.  90,  HORIZON  OBSTACLE  HEIGHT]  or 
calculated  as  indicated  in  Figure  5. 

h.  ,  Height  of  facility  horizon  [km]  above  the  reflecting  surface 

(Figure  6),  from  (39). 

h.  9  Height  [km]  of  aircraft  horizon  above  the  reflecting  surface, 

Lrv  from  (47). 

h^i  2  Horizon  elevations  [km]  above  the  msl,  from  (38,48). 

h  Elevation  [km]  of  effective  reflecting  surface  above  msl 

*  (Figure  2),  an  input  parameter  [20,  p.  83,  EFFECTIVE  REFLECTION 

SURFACE  ELEVATION]. 

* 

h  ,  p  Actual  antenna  elevations  above  reflecting  surface  elevation, 

’  from  (20),  and  shown  in  Figure  2. 

h  i  Height  [km]  of  facility  site  surface  above  msl  (Figure  2),  an 

smi  input  parameter  [20,  p.  101,  TERRAIN  ELEVATION]. 

h  ,  Height  [km]  of  facility  antenna  above  the  facility  site  surface 

(Figure  2),  an  input  parameter  [20,  p.  81]. 

h  Height  [km]  of  the  common  volume  above  the  reflection  plane 

*  for  scatter  calculations,  from  (201). 

h,  2  Antenna  elevations  [km]  above  msl,  from  (21)  and  shown  in 

'  Figure  2,  is  an  input  parameter  [20,  p.  80]. 

h^  2  Parameters  calculated  from  (80). 

H  .  Heights  [km]  defined  and  illustrated  in  Figure  12. 

C  *  C|  ,  l  >  z 

H,  9  Antenna  elevation  [km]  above  the  reflecting  surface  (Figure  8), 

from  (106). 

Hj  2  Antenna  elevations  [km]  from  (110)  and  shown  in  Figure  8. 

H^3  Significant  wave  height  [20,  p.  99]  used  in  (125). 

H  i  2  Heights  [km]  used  in  Figure  1 2  and  defined  for  different  path 

^  ’  types  in  Section  9. 

i  An  index  for  specific  transmission  loss  levels  used  in  the  dis¬ 

tribution  mixing  process  (Section  10.6).  It  has  values  from 
1  to  M. 
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IT S - FAA - 1 973  propagation  model. 

JTS-FAA-1 977  propagation  model. 

Institute  for  Telecommunication  Sciences. 

/-T  as  in  (154)  or  index  (1  to  N)  for  specific  transmission 
loss  distributions  used  in  the  distribution  mixing  process 
(Sec.  10.6). 

An  adjusted  earth  radius  factor,  from  (103). 

The  ratio  [dB]  between  the  steady  component  of  received  power 
and  the  Rayleigh  fading  component  that  is  used  to  determine 
the  appropriate  Nagagami-Rice  distribution  [34,  p.  V-8]  for 
Y7T(q),  from  (258). 

Rounded  earth  diffraction  parameters,  from  (65). 

Rounded  earth  diffraction  parameters,  from  (67). 

Kt  values  in  the  1 ine-of-sight  region,  from  (254). 

K  va’ues  at  the  radio  horizon;  i.e.,  K.  nc  at  d  =  dul  , 
from  (254).  L0S  ML 

K  value  associated  with  tropospheric  multipath,  from  (256). 


^234  Rounded  earth  diffraction  parameters,  from  (66). 

Total  ray  length  [km]  in  scatter  calculation,  from  (204). 
og  f jmmon  (base  10)  logarithm. 

oge  Natural  (base  e)  logarithm. 


Ray  lengths  [km]  to  the  cross-over  point  of  the  common  volume 
in  tropospheric  scatter  calculations,  from  (203). 


OS  Une-of- sjght. 

(q)  Transmission  loss  [dB]  values  not  exceeded  during  a  fraction 

q  of  the  time.  These  values  may  represent  instantaneous  levels 
depending  upon  the  time  availability  option  selected  (Sec.  10), 
and  are  calculated  using  (1). 

Basic  transmission  loss  [dB]  for  the  free  space,  from  (226). 


76 


br 

Lb(0.5) 

L9P 

mhos 

msl 

M 

MHz 

Md 


M 


Ka 


M. 


Mr 


M. 


n 


1,2,3 


N 

N 

N-uni ts 


A  reference  level  of  basic  transmission  loss  [dB],  from  (3). 
Median  basic  transmission  loss  [dB],  from  (2). 

Loss  [dB]  in  path  antenna  gain  used  in  and  discussed  after  (3); 

current  model  assumes  L  „  =  0. 

9P 

Unit  of  conductance  or  siemens. 

Mean  sea  level . 

Number  of  transmission  loss  levels  used  in  mixing  distributions 
which  is  also  the  final  value  for  the  index  i  (Sec.  10.6). 

Megahertz  (10^  Hz). 

Slope  [dB/km]  of  combined  diffraction  attenuation  line  for 
beyond- the- horizon,  from  (98). 

Slope  [dB/km]  of  rounded  earth  diffraction  attenuation  line  for 
path  F-G-ML  (Figure  7),  from  (55,77). 

Slope  [dB/km]  of  the  K  value  line  used  just  beyond  the  radio 
horizon,  from  (255). 

Slope  [dB/km]  of  the  diffraction  attenuation  line  used  just 
inside  the  radio  horizon,  from  (176). 

Slope  [dB/km]  of  rounded  earth  diffraction  line  attenuation  for 
path  0-A  (Figure  7),  from  (55,77). 

Slope  [dB/km]  of  rounded  earth  diffraction  attenuation  line  for 
path  p,  from  (55,77). 

Slope  [dB/km]  of  successive  scatter  attenuation  points  (A  versus 
d  line),  discussed  after  (221). 

A  power  used  in  the  ionospheric  scintillation  frequency  scaling 
factor  (Sec.  10.5),  from  (262). 

Parameters  for  (234,235)  that  are  discussed  following  (2351. 

Refractivlty  [N-units]  for  a  height  h  in  an  exponential  atmos 
phere,  from  (17). 

Number  of  distributions  to  be  mixed  which  is  also  the  final 
value  for  the  index  j  (Sec.  10.6). 

Units  of  refractivity  [3,  Sec.  1.3]  correspondi ng  to  10^  (re¬ 
fractive  index  -1 ) . 
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Pa(q) 


qll,ij,MN 


Qf\2 


Qol,2 


Rational  Telecommunications  and  information  Administration. 

Minimum  monthly  mean  surface  refractivity  (N-units)  referred  to 
.iisl,  an  input  parameter  [20,  p.  94]. 

Minimum  monthly  mean  surface  refractive  (N-units)  at  effective 
reflection  surface,  calculated  from  NQ  via  (14). 

Power  [dBW]  that  is  available  at  the  terminal  of  an  ideal 
(loss-less)  receiving  antenna  for  at  least  a  fraction  q  of 
the  time,  from  (5). 

Total  power  [dBW]  radiated  from  the  transmitting  antenna,  used 
in  (7). 

Dimensionless  fraction  of  tinje  used  in  time  availability  speci¬ 
fication;  e.g.,  in  L.(O.l),  q  =  0.1  implies  a  time  availability 
of  10  percent.  D 

Parameters  used  in  tropospheric  scatter  calculations,  from  (216) 

Time  availabilities  for  mixed  distributions  that  correspond  to 
specific  transmission  loss  levels,  from  (264). 

rune  availabilities  for  each  transmission  loss  level  (index  i) 
of  each  transmission  loss  distribution  (index  j)  involved  in 
the  distribution  mixing  distribution  (Sec.  10. S). 

Parameters  used  in  tropospheric  scatter  calculations,  from  (189) 

Parameters  used  in  tropospheric  scatter  calculations,  from  (193) 

Parameters  used  in  tropospheric  scatter  calculations,  from  (190) 

Parameters  used  in  tropospheric  scatter  calculations,  from  (194) 

Parameters  used  in  tropospheric  scatter  calculations,  from  (188) 

Ray  length  [km]  used  in  the  calculation  of  free  space  loss, 
from  (225). 

Ray  length  [km]  for  beyond- the-horizon  paths,  from  (224). 


r£  A  distance  [km],  from  (131). 


r  Effective  i  ay  lengths  [km]  for  attenuation  associated  with 

’  ’  oxygen  absorption,  (r  )  rain  storm  attenuation  (rgs),  and 

water  vapor  absorption  (r  ),  from  (227). 

cW 
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r,  i  0  Antenna  to  horizon  ray  lengths  [km]  for  airless  earth,  from 

Ll,d  (223). 

rQ  Direct  ray  length  [km]  shown  in  Figure  8,  from  (112). 

r  In-storm  ray  length  [km]  used  in  rain  attenuation  calculation, 

b  from  (259). 

rWH  Within-the-horizon  ray  length  [km]  for  airless  earth,  from  (222). 

r,  ?  Segments  of  reflected  ray  path  shown  in  Figure  8  and  components 

*  of  r^.  ^rom  (^0). 

r-j ^  Reflected  ray  path  length  [km]  as  shown  in  Figure  8,  from  (113). 

rleos  w  Partial  effective  ray  lengths  [km]  for  oxygen  layer,  rain  storm 

'  j  or  water  vapor  layer  or  calculated  using  the  relationships 

given  in  Figure  12,  as  indicated  in  Section  9. 

r2eo,w 

R  Magnitude  of  complex  plane  earth  reflection  coefficient,  from 

(155). 

RAYTRAC  A  computer  subroutine  used  to  trace  a  ray  through  an  exponential 

atmosphere  [14,  p.  182]. 

Rc  Magnitude  of  effective  reflection  coefficient  for  counterpoise 

,y  or  ground  reflection,  from  (156). 

R.  Diffuse  component  of  surface  reflection  multipath,  from  (129, 

251). 

Rg  |  Magnitudes  of  complex  plane  earth  reflection  coefficients  for 

’  ’  eliptical,  horizontal,  and  vertical  polarization. 

R  .  Partial  effective  ray  lengths  [km]  inside  oxygen  layer,  rain- 

storm  or  water  vapor  layer,  as  calculated  as  shown  in  Figure  12. 

R  A  parameter  used  in  the  calculation  of  the  divergence  factor, 

from  (121 ) . 

Rs  Specular  component  of  surface  reflection  multipath,  from  (250). 

RT  Magnitude  of  adjusted  effective  reflection  coefficient  for 

y,c  earth  or  counterpoise  reflection,  from  (119). 

s  Modules  of  asymmetry  used  in  tropospheric  scatter  calculations, 

from  (209) 
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SHF 


Super  High  frequency  [3  to  30  GHz]. 

SirT^  Inverse  sine  with  principal  value. 

S  Scattering  efficiency  term  [dB]  used  in  tropospheric  scatter 

calculations,  from  (208). 


SR(q)  Power  density  at  receiving  antenna  locations  [dB- W/sq  m]  for 

at  least  a  fraction  q  of  the  time,  from  (8). 
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A  Fresnel  integral  [34,  p.  1 1 1-18] •  for  (92,137,167). 

Scattering  volume  term  [dB]  of  tropospheric  scatter  calculations, 
from  (219). 

A  parameter  defined  and  used  in  (83),  the  G_  formulation. 

h 

Relaxation  time  [ps]  used  in  the  calculation  of  water  surface 
constants,  for  (151). 

Inverse  tangent  [rad]  with  principal  value. 

A  program  name  [20,  p.  9]. 

Storm  height  or  layer  thickness  (Figure  12)  jsed  in  attenuation 
calculations  for  oxygen  absorption  (Tg0),  rain  storm  attenuation 

(T  ),  or  water  vapor  absorption  (r). 

GS  GW 

Ultra-figh  frequency  [300  to  3000  MHz], 

Knife-edge  diffraction  parameter,  from  (90). 

Knife-edge  diffraction  parameters  used  to  determine  f  c* 
from  (135,165).  9*  * 

Very  High  frequency  [30  to  300  MHz). 

Volts  per  volt. 

A  parameter  [dB],  from  (235,246). 

Variability  for  specific  climate  or  time  block,  from  (265). 

Variability  adjustment  term  [dB],  from  (240). 

Variability  levels  (V, ,  ...  V-,  ...  VM)  used  in  mixing  distri¬ 
butions  (Sec.  10.6).  ‘ 

A  weighting  factor  used  in  combining  knife-edge  and  rounded 
earth  diffraction  attenuations,  from  (94). 


W  Watts. 


'R 

fR0 

'l,2 
1  >  j » N 

a 

A1 ,2 
vl  ,2 

1,2, 3,4 

1,2 

(q) 

(0.1) 

or 

(0.9) 

c 

DU^q) 

e(q) 

«i(q) 


A  relative  power  level  for  the  Rayleigh  fading  component 
associated  with  tropospheric  multipath,  from  (257). 

A  relative  power  level  for  the  Rayleigh  fading  component 
associated  with  surface  reflection  multipath,  from  (252). 

A  relative  power  level  associated  with  the  ray  optics  formula¬ 
tion  used  in  the  line-of-sight  region,  from  (174). 

Parameters,  from  (72). 

Weighting  factors  (W, ,  ...  W.,  ...  WM)  used  in  mixing  distri¬ 
butions  (Sec.  10.6).  J  M 

A  parameter,  from  (64). 

A  parameter  used  in  tropospheric  scatter  calculations,  from 
(198). 

Parameters  used  in  tropospheric  scatter  calculations,  from  (202) 

Parameters  used  in  tropospheric  scatter  calculations,  from 
(212,213). 

Parameters  [km],  from  (69,70). 

Parameters  [dB],  from  (73). 

Variability  [dB  greater  than  median]  of  hourly  median  received 
power  about  its  median,  from  (237). 

A  reference  variability  level  of  hourly  median  received  power 
about  its  median  and  is  used  to  calculate  the  other  percentiles, 
from  (236). 


A  complex  parameter  used  in  the  calculation  of  the  plane  earth 
reflection  coefficient,  from  (153). 

Total  variability  [dB  greater  than  median]  of  D/U  about  its 
median,  Ypy(0.5)  =  0,  where  q  is  the  fraction  of  time  for  which 

a  particular  value  is  exceeded. 

Effective  variability  [dB  greater  than  median]  of  hourly  median 
received  power  about  its  median,  from  (244,245). 

Initial  value  of  Y  (q),  from  (239). 


Vj(q) 

v,c(q) 

Vr(q) 

v0(°.D 

or 

Yo(0.9) 

YT 

Yv 

Y136(q) 

Y,(q) 

Y,(q) 


zal  ,2 
2bl  ,2 
zl,2 
Zal  ,2 
ZA1 , 2 
Zbl ,  2 


Variability  [dB]  associated  with  ionospheric  scintillation, 
from  (263). 

Y I (q ) | c  for  a  particular  distribution  to  be  used  in  the  mixing 
process  to  obtain  resultant  Yj(q),  from  (267). 

Variability  [dB  greater  than  median]  associated  with  rain 
attenuation,  from  (261). 

A  reference  variability  level  used  to  calculate  Y(0.1)  or  Y(0.9), 
from  (235,246). 


A  top  limiting  level  when  the  calculations  are  in  the  lobing 
mode  in  the  1 ine-of-sight  region,  from  (241). 

A  parameter  from  (133). 


Variability  [dB]  for  (263)  that  is  associated  with  ionospheric 
scintillation  at  136  MHz  [15,  p.  45], 

Variability  [dB  greater  than  median]  of  received  power  used  to 
describe  short-term  (within-the-hour)  fading  associated  with 
multipath  where  q  is  the  faction  of  time  during  which  a  particu¬ 
lar  level  is  exceeded  (Sec.  10.3). 

Total  variability  [d3  greater  than  median]  of  received  power 
about  its  median,  Y£(0.5)  =  0,  where  q  is  the  fraction  of  time 
for  which  a  particular  value  is  exceeded,  from  (229).  These 
values  may  represent  instantaneous  levels  or  hourly  median  levels 
depending  upon  the  time  availability  option  selected  [20,  p.  103, 
TIME  AVAILABILITY  OPTIONS]. 


A  parameter  from  (102). 

Parameters  used  in  tropospheric  scatter  calculations,  from  (186). 


Parameters  used  in  tropospheric  scatter  calculations,  from  (187). 
Parameters  [km],  from  (107). 

Parameters  used  in  tropospheric  scatter  calculations,  from  (191). 


Parameters  used  in  tropospheric  scatter  calculations,  from  (195). 
Parameters  used  in  tropospheric  scatter  calculations,  from  (192). 


An  angle  [rad]  shown  in  Figure  8,  from  (111). 

An  angle  [rad]  used  in  Figure  12  and  defined  for  different 
path  types  in  Section  9. 

A  parameter  [per  km]  used  in  tropospheric  scatter  calculations, 
from  (207). 

A  parameter  in  tropospheric  scatter  calculations,  from  (185). 

Surface  absorption  rates  [dB/km]  for  oxygen  or  water  vapor  for 
(228). 

A  parameter  used  in  the  calculation  of  the  surface  roughness 
factors,  from  (126). 

Terrain  parameter  [km]  used  to  characterize  terrain  irregularity 
It  is  a  model  input  parameter  [20,  p.  101,  TERRAIN  PARAMETER]. 

Adjusted  effective  altitude  correction  factors,  from  (105). 

Interdecile  range  of  terrain  heights  [m]  above  and  below  a 
straight  line  fitted  to  elevations  above  msl ;  estimated  from 
(124)  which  is  based  on  previous  work  [26,  eq.  3]. 

Effective  altitude  correction  factors,  from  (25). 

Ah  expressed  in  meters  for  (124). 

Refractivity  gradient  [N-units/km]  used  in  defining  exponential 
atmospheres,  from  (19). 

Path  length  difference  [km]  for  rays  shown  in  Figure  8 
( r i 2  -  rQ),  from  (101,114). 

/ 

Ar  [km]  for  earth  or  counterpoise  reflection,  from  (114). 

Dielectric  constant,  model  input  parameter  [20,  p.  99,  SURFACE 
TYPE  OPTIONS]. 

Complex  dielectric  constant,  from  (152). 

Dielectric  constant  representing  the  sum  of  electronic  and 
atomic  polarizations.  For  water,  c  =  4.9  in  (150,151). 

Static  dielectric  constant  for  water  [15,  p.  26]. 


Parameters  used  in  tropospheric  scatter  calculations,  from 
(205,206). 
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A  parameter  used  in  tropospheric  scatter  calculations,  from 

(211). 

Scattering  angle  [rad]  used  in  tropospheric  scatter  calculations. 
It  is  the  angle  between  transmitter  horizon  to  common  volume  ray 
and  the  common  volume  to  receiver  horizon  ray  as  both  leave 
their  crossover  point,  from  (197). 

Angles  [rad],  from  (178). 


Angles  [rad]  used  in  tropospheric  scatter  calculations;  angles 
between  the  common  volume  horizontal  and  the  horizon  to  the 
common  volume  rays  as  they  leave  their  crossover  point, 
from  (196). 

An  angle  [rad]  shown  in  Figure  9,  from  (130). 

An  initial  value  for  the  facility  horizon  elevation  angle  [rad] 
that  is  based  on  effective  earth  radius  geometry,  and  is  shown 
in  Figure  4.  It  may  be  specified  [20,  p.  90,  HORIZON  OBSTACLE 
ABOVE  HORIZONTAL  AT  FACILITY]  or  calculated  as  indicated  in 
Figure  5. 

Smooth  earth  horizon  ray  elevation  angle  [rad]  at  the  appropriate 
terminal  as  determined  from  ray  tracing  (Sec.  3.1).  Illustrated 
in  Figure  3. 

Final  value  for  smooth  earth  horizon  ray  elevation  angle  [rad] 
at  the  appropriate  terminal,  from  (30). 

Elevation  angle  [rad]  of  horizon  from  the  facility,  from  (40). 

Elevation  angle  [rad]  of  aircraft  horizon  ray  snown  in 
Figure  6,  from  (49). 

An  angle  [rad],  from  (163). 

Magnitude  of  earth  facility  latitude  [deg]  for  (262). 

Elevation  angles  [rad]  of  the  ground  reflected  rays  at  the 
terminal  antennas,  from  (148). 

Parameter  [rad]  used  to  calculate  6^  ^  from  (115). 


An  angle  [rad],  from  (158). 

Direct  ray  elevation  angles  [rad]  at  the  terminal  antennas, 
from  (147). 


84 


-  ' l  Si  .  "  -  *'■  ■ 


Angles  [rad]  used  in  (135),  from  (132,134). 

Elevation  angle  [rad]  of  horizon-to-aircraft  ray  at  facility 
horizon  (Figure  4),  from  (42). 

Horizon  elevation  angle  [rad]  adjustment  terms,  from  (146). 
Central  angle  [rad],  from  (117). 

Angles  [rad]  used  in  (148),  from  (116). 

An  angle  [rad],  from  (180). 

Smooth  earth  horizon  ray  elevation  angles  [rad]  that  a^e 
obtained  using  ray  tracing  horizon  distances  with  an  effective 
earth  formulation  and  shown  in  Figure  3,  from  (23). 

Central  angles  [rad]  below  the  smooth  earth  terminal  to 
horizon  distances,  from  (29). 

Diffraction  angle  [rad]  for  the  0-A  path  (Figure  7),  from  (89). 

Angles  [rad]  shown  in  Figure  8,  from  (108). 

Angles  [rad],  from  (57,58). 

An  angle  [rad],  from  (159). 

An  angle  [rad],  from  (164). 

Wave  number  [per  km^,  from  (214). 

Wavelength  [km],  from  (10). 

Wavelength  [m'j,  from  (10). 

Microseconds  [10‘6  sec]. 

The  constant  3 . 141 592654 . 

Parameters  used  in  tropospheric  scatter  calculations,  from 
(215). 

Surface  conductivity  [mho/m],  a  model  input  parameter  [20,  p.  99, 
SURFACE  TYPE  OPTIONS]. 


The  root-mean- square  deviation  [m]  of  terrain  and  terrain 
clutter  within  the  limits  of  the  first  Fresnel  zone  in  the 
dominant  reflecting  plane;  estimated  from  (125),  which  is 
based  on  previous  work  [26,  Annex  3,  eqs.  3.6a,  3.6b]. 

Ionic  conductivity  [mho/m]  for  water,  [15,  p.  26]. 

Phase  advance  associated  with  complex  earth  reflection  coeffi¬ 
cient  used  in  (155). 

Phase  shift  [rad]  of  effective  reflection  coefficient  for  earth 
or  counterpoise  reflection,  used  in  (156). 

Knife-edge  diffraction  phase  shift  [rad]  for  earth  or  counter¬ 
poise  reflection,  from  (138). 

Total  phase  shift  [rad]  of  effective  reflection  coefficient  for 
earth  or  counterpoise  reflection,  from  (157). 

Grazing  angle  shown  in  Figures  8  and  9. 

F.xpression  evaluated  for  specific  conditions  such  as  climate 
or  time  block  in  (267). 
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